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1  1 

SUMMARY 


This  report  provides  e  deteiled  description  of  the  methods 
used  in  the  computer  program  CHAOS  for  calculating  radar  cross-sections. 


1 .  INTRODUCTION 


The  object  of  this  report  is  to  provide  a  detailed  description 
of  the  methods  used  in  the  CHAOS  computer  program  for  calculating  radar 
cross-sections.  The  program  is  written  in  FORTRAN. 


The  various  physical  approximations  aTe  described,  4» 
section  -3^  These  help  to  simplify  Maxwell's  equations.  At  first  sight 
these  would  appear  to  limit  the  scattering  bodies  to  open  wire 
structures;  however,  extensions  to  surfaces  both  open  and  closed  follow 
by  further  approximation  (section  S>.  The  numerical  methods  used  to 
solve  these  equations  are  described  in  oeetien  3. 

'r\ 

The  problem  of  wire  junctions  is  tackled  in  section  4  while 
section  6  deals  with  finite  conductivity  and  shows  how  lumped  loads  may 
be  attached  to  the  structure.  Some  special  techniques  for  dealing  with 
bodies  having  rotational  or  left-right  symmetries  are  described  in 
section  7.  Sections  8  and  9  briefly  describe  the  sparse  matrix 
approximation  and  the  complex  frequency  version  of  CHAOS  respectively. 

Three  appendices  are  also  provided.  Appendix  A  deals  with  the 
formulae  used  <$o  calculate  the  impedance  matrix  elements  and  gives 
reasons  for  various  restrictions  regarding  segment  length  and  radius. 
Appendix  B  describes  the  conventions  adopted  throughout  the  programme 
concerning  coordinate  systems  and  polarisation.  A  list  of  frequently 
used  symbols  is  given  in  appendix  C. 


2. 


THE  DERIVATION  OF  THE  CHAOS  EQUATIONS 


2.1  Maxwell *8  equations  in  the  frequency  domain 

We  start  from  Maxwell's  equations.  These  are:- 
Faraday's  law  of  induction: 


Gauss'  law: 

V.D  *  0, 

Maxwell's  generalisation  of  Ampere's  law: 
3D 

^  “  ST  +  * 

and  finally 


....(1) 


(2) 


.(3) 


V.B  -  0. 

In  these  equations  the  symbols  have  the  following  meaning:- 


(4) 


E  *  complex  electric  field, 

D  *  complex  electric  flux, 

H  *  complex  magnetic  field, 

B  *  complex  magnetic  flux, 

J  *  complex  conduction  current, 
p  *  complex  electric  charge, 
t  *  time, 


and  the  rationalised  system  of  MKS  units  has  been  used.  At  this  point  we 
have  the  choice  of  working  in  the  time  domain,  so  that  we  can  deal  with  an 
initial  pulse  which  is  some  arbitrary  function  of  time,  or  of  specifying 
some  particular  time  behaviour.  We  choose  the  latter  and  assume  a  time 
variation  of  eia)t  where  u)  is  the  frequency  of  the  electromagnetic  wave. 
This  is  known  as  working  in  the  frequency  domain.  Thus,  this  type  of 
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approach  is  ideally  suited  to  calculating  radar  cross-sections  for 
radars  of  the  continuous  wave  type.  For  radars  which  send  out  short 
pulses  containing  components  spread  over  many  frequencies ,  it  is 
necessary  to  Fourier  analyse  the  pulse  and  do  calculations  for  all  the 
frequencies  in  the  pulse  before  combining  together  to  get  the  scattered 
field  and  radar  cross-section.  Despite  the  apparent  clumsiness  of  this 
approach,  it  appears  to  be  more  reliable  than  working  directly  in  the 
time  domain. 

By  specifying  that  all  the  quantities  in  Maxwell's  equations 
have  a  time  variation  e  ,  we  obtain 

-  id)B  } 

"  } 

P  } 

}  ....(5) 

iwD  +  J  } 

"  “  ) 

0.  } 

2.2  Homogeneous  materials  assumption 

There  is  a  relation  between  the  variables  D  and  E  and  another 
relation  between  the  variables  B  and  H.  It  is  known  that  for  some 
materials  these  pairs  of  variables  are  linearly  related  and  we  have 

D  =  GE, 

B  =  MH, 

where  e  and  y  are  scalars  called  the  dielectric  constant  and  perme¬ 
ability  respectively.  For  homogeneous  materials  e  and  y  are  independent 
of  position  in  the  body.  So  equations  (5)  reduce  to 

7^E  *  -  iwyH  } 

} 

V.J  *  p/e  }  ....<« 

7/sH  »  iweE  +  J  } 

"  } 

7.H  -  0.  } 


7/.E  = 
7.D  - 
7~H  = 
V.B  = 
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2.3  Helmholtz 1 s  equations  for  potentials 

From  the  last  of  equations  (6)  we  have 

H  =  jj  V.A,  ....(7) 

where  A  is  the  magnetic  vector  potential.  The  first  of  equations  (6)  may 
now  be  integrated  to  give 

E  =  -  iwA  —  V4> ,  ....  (8) 

where  $  is  the  electric  scalar  potential  and  is,  at  present,  some 
arbitrary  function.  By  substituting  in  the  remaining  two  equations  we  get 
the  Helmholtz  equations 


V2A  +  k2  A  =  -  yj, 


V2<J)  +  k2<J>  =  -  p/e, 


(9) 


where  k2  =  u)2ey  and  we  have  chosen  to  impose  on  $  the  so  called  Lorentz 
condition 


V.A  *  -  iwEy<J).  ....  (10) 

2. A  Solution  of  Helmholtz's  equations  for  conductors 

Now  we  can  formally  solve  the  Helmholtz  equations  and  since  for 
conductors  it  is  known  that  the  current  J  and  charge  p  reside  on  or  very 
close  to  the  surface,  the  vector  equation  gives 

■~ik.R 

A(r)  =  y  #  J(r')  ds’,  ....(11) 

where  R  is  the  distance  between  the  observation  point  r  and  some  point  r' 
on  the  surface  of  the  conductor.  The  scalar  equation  gives 

1  r/  -ikR 

-Mr)  -  -  &  P(r')  £_  ds'.  ....(12) 
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In  both  these  expressions  ds1  is  an  element  of  area  of  the  surface  of  the 
conductor  and  the  integrations  are  over  the  whole  surface.  There  is  one 
further  equation  relating  the  current  J  and  the  charge  P  which  is  the  law 
of  conservation  of  charge  obtained  by  differentiating  the  third  of 
equations  (5) 


V.J  3  -  icop. 


....(13) 


2.5  Boundary  conditions 

At  a  boundary  between  two  different  materials,  the  following 
conditions  are  observed 


(E(l)  -  E(z))^n  -  0, 
n.(D(l)  -  D(2J)  =  p 

iu(H(l)  -  H(2))  =  J  , 
—  —  —  — s 

n.(B(l)  -  B(2))  =0, 


....(14) 


where  the  superscripts  (1)  and  (2)  refer  to  the  fields  on  either  side  of 
the  material  boundary,  and  n  is  the  unit  normal  to  the  surface  pointing 
into  region  (1).  Here  pg  and  Jg  are  electric  surface  charge  and  current 
respectively.  Since  we  are  solving  for  the  electric  field  it  is  the  first 
two  which  concern  us. 


2.6  Thin  wire  approximation 

The  integrals  in  (11)  and  (12)  can  be  further  simplified  by 
assuming  that  we  are  dealing  with  thin  wires  rather  than  surfaces.  While 
this  may  appear  to  restrict  the  program  to  dealing  with  wire  structures 
only,  we  shall  see  in  section  5  that  surfaces  can  be  modelled  by  wire 
meshes  and  this  appears  to  be  the  best  method  available  for  dealing  with 
arbitrarily  shaped  surfaces  of  finite  conductivity. 
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Following  from  the  thin  wire  assumption,  we  assume  that  the 
circumferential  current  is  negligible  and  replace  the  axial  surface 
current  and  the  surface  charge  by  a  line  current,  I,  flowing  along  the 
wire  axis  and  a  line  charge  o.  So  we  have 


I_  *  2iraJ, 
o  =  2nap, 


where  a  is  the  wire  radius.  This  leads  to 

-ikR 

A(r)  =  y  /  1(4 ' )  —u  d4  ’ , 


wire 


1  ,  -ikR 

<J.(r)=i  /  a(l')  Vrr^', 

e  wire 


o(4)  =  -  V.I(4)  =  i  1(4), 
go  -  to  d4  - 


....(15) 


(16) 


....(17) 


....(18) 


,(19) 


where  4  is  the  measure  of  length  along  the  wire.  We  only  try  to  satisfy 
the  boundary  condition  on  the  axial  component  of  the  electric  field:- 


(E^l}  -  E^z))  =  0, 
-tan  -tan  ’ 


•  •  •  • 


(20) 


where  the  subscript  tan  denotes  the  component  of  the  electric  field 
parallel  to  the  wire  axis. 


2.7  Currents  due  to  plane  wave  excitation 

Let  us  now  assume  that  the  scattering  body  will  be  illuminated 
by  a  plane  electric  wave,  E1,  of  known  amplitude  and  polarisation.  Then 
the  total  electric  field  may  be  written 

E  -  E1  +  ES,  ....(21) 

s 

where  E  is  the  unknown  electric  field  scattered  from  the  body.  Bearing 
in  mind  equations  (8),  (17)  and  (18)  let  us  write 
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E8  *  -  iu)A8  -  7<J>9 , 


f 

J 

-ikR 

AS(r)  -  y 

I(i’) 

e 

4itr 

dl 

wire 

4>S(r)  =  — 
Y  -  e 

-ikR 

/ 

a(V) 

e 

dQ 

wire 

4ttr 

....(22) 

....(23) 

....(24) 


It  will  be  noticed  that  the  superscript  s  has  been  omitted  from  I  and  o  in 
equations  (23)  and  (24).  This  is  because  a  plane  wave  satisfies  Maxwell's 
equations  (6)  with  p  =  0  and  J  =  0,  so  that  the  whole  of  the  current  and 
charge  distributions  induced  on  the  body's  surface  by  the  incident  plane 
wave  contribute  to  the  scattered  fields  as  one  would  expect. 


The  boundary  condition  (20)  may  be  written 


<lia»  *  -  — tan^  *  °- 


....(25) 


2.8 


Equation  for  the  currents 


We  use  equation  (25)  because  this  relates  the  scattered  field 

on  the  body  surface  to  the  incident  field  E1  and  the  field  inside  the 

(2)  (2)  “ 

wire  E  .  For  a  perfect  conductor  E  is  zero  while  for  a  wire  of 

finite  conductivity,  it  may  be  calculated  approximately  and  in  section  5 

an  expression  of  the  form 


I 

u>- 


(26) 


will  be  derived  where  is  the  known  impedance  matrix  of  the  wires. 


So  by  substituting  equations  (19),  (22),  (23),  (24)  and  (26) 
into  equation  (25)  we  obtain  an  equation  for  the  currents  I  of  the  form 

.did') 

%  I  _ 

....(27) 


Z  I  -  E. 
u>—  —tan 


-  i  / 

wires 


1  M 


-ikR 
e _ 

4ttR 


dr 


This  equation  is  sometimes  referred  to  as  the  potential  integro- 
differential  equation  to  distinguish  it  from  other  variants  such  as 
Pocklington's  equation,  the  magnetic  vector  potential  equation  and 
Ha lien's  equation. 
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3. 


NUMERICAL  REPRESENTATION  OF  THE  PHYSICS 


3. 1  Representation  of  the  currents 

We  now  expand  the  currents  in  terms  of  a  finite  series  of  known 
basis  functions  with  unknown  amplitudes.  If  we  were  dealing  with  a  single 
wire,  the  obvious  choice  would  be  a  set  of  whole  wire  functions  like,  for 
instance,  a  Fourier  series  expansion.  However,  we  wish  to  be  able  to  deal 
with  bent  wires  and  arbitrarily  shaped  wire  meshes  so  it  is  more  useful 
to  divide  the  wires  up  into  a  series  of  short  segments  and  define  a  set  of 
functions  each  of  which  is  non-zero  over  a  few  segments  only. 


We  choose  triangle  functions  as  used  by  Chao  and  Strait  [1], 
Mautz  and  Harrington  [2]  and  Kuo  and  Strait  f3],  As  shown  in  figure  1(a), 
each  triangle  function,  T^,  is  based  on  four  segments  and  overlaps  its 
neighbour  by  two  segments.  It  is  of  unit  height.  Thus,  we  expand  the 
current  as 

N 

I(Jl)  =  y  I.T.U),  - (28) 

i-i  1_1 

where  N  is  the  number  of  triangle  functions  needed  to  cover  the  wires.  In 
this  way  a  piecewise  linear  approximation  to  the  currents  is  obtained  as 
depicted  in  figure  1(b).  The  amplitudes  1^  are  the  unknowns  which  we 
require  to  find. 

Other  commonly  used  basis  functions  are  step  functions, 
sinusoidal  functions  and  a  three  term  function  consisting  of  constant, 
sine  and  cosine  terms. 


3.2  Choice  of  weighting  function 

Since  we  have  chosen  a  finite  number  of  basis  functions  T.  (i  ■ 
1,  ...  N),  we  need  to  derive  N  equations  from  equation  (27)  to  enable  us 
to  solve  for  the  unknown  amplitudes  1^.  We  do  this  by  multiplying 
equation  (27)  in  turn  by  N  weighting  functions  and  integrating  over  the 
wire  length.  These  weighting  functions  could  be  the  delta  function 
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positioned  at  the  centre  of  each  triangle  function  and  in  this  way  we 
would  satisfy  equation  (27)  exactly  at  these  points.  However,  the 
solution  could  be  badly  in  error  away  from  these  points  and  it  has  been 
found  better  to  use  the  basis  functions  as  weighting  functions.  In  this 
way  we  do  not  satisfy  equation  (27)  exactly  anywhere,  except  by  accident, 
but  we  minimize  the  overall  error  of  the  solution  for  the  current  in  the 
least  squares  sense.  This  latter  scheme  is  known  as  the  Galerkin 
procedure,  while  the  former  is  called  point  matching  or  collocation. 
There  is  an  added  attraction  to  the  Galerkin  procedure  in  that  the 
resulting  matrix  of  impedances  is  symmetric  as  will  be  seen  shortly. 
Thus,  only  half  the  matrix  need  be  stored  in  the  computer  and  the 
principle  of  reciprocity  is  automatically  satisfied.  This  is  not  the 
case  for  any  other  set  of  weighting  functions. 


3.3  Form  of  impedance  matrix 


By  substituting  (28)  into  (27)  and  carrying  out  the  Galerkin 
procedure  we  arrive  at  the  set  of  equations 


N 


y  Z..I.  =  V.,  j  -  1  ...  N, 

.L  Hi  1  •*  1 

1=1  J  J 

where  the  impedance  matrix  is  given  by 

Zii  “  {  <&Z  .T.(4).Ti(Jt)  +  ic»M  /  dl  f  dZ' T.U).T.(i') 

wires  J  wires  wires  “-1  ‘“1 


....(29) 


n  e 


*Se  (  jV*'1  e2iR 

Wires  J  ui  ra a  *  4lTR 


-ikR 

4ttR 


wires 


wires 


....(30) 


and  the  excitation  voltages  are  given  by 


V  -  J  (D.jj 

J  wires  J 


....(31) 


From  the  form  of  the  last  term  in  equation  (30)  it  appears  that  Z..  is  not 

ji 

symmetric.  However,  we  can  manipulate  the  last  term  to  obtain 
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dT.(i')  -ikR 

k  !  dlV°-7  {  il'  -Jr-  ' 

Wires  J  «1  ret  a 


wires 


4*R 


i  r  dr  ".(£')  -ikR 

Sc  I  “V11!  dt'-3^-W- 

Wires  J  vi  ros  ‘♦7TK 


vires 


dT.(t’)  ,-ikR|l-l„  ;  H  ...  .•  ...  "‘i™  "  «  ~ 

1-  fT. a)  /  <*•  -557“^  _  -  k  f.  dt  Si  Tia)  !  “  ~nr~  Ta¬ 
me  |  j  wires  41,11 1 £-4  *•  wires  wires 


dT.(t’)  -i»* 


Now  the  first  term  is  zero  because  the  currents  at  the  free  ends  of  the 

wires  are  assumed  to  be  zero.  Thus,  Z..  can  be  written 

*  Ji 


Z..  =  /  d£Z  .T.a).T.(jl) 


n 


wires 


IDl-J  “V 


.if  Mf  dfLlT  (D.T^l’)  -  T«)  Ijt  T.(f)| 

wires  wires  J  J  J 


....(32) 
-ikR 


4irR  ’ 


which  is  symmetric  in  i  and  j  which  means  that  the  principle  of 
reciprocity  is  automatically  satisfied.  At  first  sight  it  may  appear  that 
equation  (32)  is  only  valid  for  a  set  of  unconnected  wires  with  both  ends 
free.  However,  once  the  method  of  joining  wires  together  is  described  (in 
section  4)  it  will  become  clear  that  (32)  holds  for  any  wire  structure. 


It  is  perhaps  worth  pointing  out  that  if  we  were  using  a 
surface  patch  approximation  rather  than  the  thin  wire  approximation,  the 
impedance  matrix  analogous  to  that  given  by  equation  (30)  would  be 
unsymmetric  even  if  Galerkin's  technique  were  used.  It  is  not  possible  to 
devise  an  equation  analogous  to  equation  (32)  because  the  current  at  the 
edge  of  a  patch  is  not  necessarily  zero.  Thus,  to  get  a  symmetric  matrix 
one  must  work  this  in  from  the  start  by  using  the  principle  of 
reciprocity.  This  is  the  approach  adopted  by  Richmond  [4]  in  his  thin 
wire  treatment. 


In  CHAOS  the  voltages  given  by  equation  (31)  are  for  plane  wave 
excitation.  It  would  not  be  difficult  to  modify  this  equation  to  cope 
with  a  point  source,  as  is  seen  in  antenna  problems,  so  that  radiation 
patterns  could  be  calculated  rather  than  radar  cross-sections.  However, 
this  has  not  been  done. 
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3.4  Matrix  factorization 

Equation  (29)  £or  the  unknown  current  amplitudes  is  solved 
using  the  Harwell  matrix  factorisation  routine  MA29  [5]  modified  for 
complex  elements.  The  method  is  due  to  Fletcher  [6)  and  the  factorisation 
takes  the  form 

T 

Z  -  LDL  ,  ....(33) 

where  L  is  lower  triangular  and  D  is  block  diagonal,  the  blocks  not 
exceeding  2  x  2  in  size.  As  a  check  on  the  reliability  of  the  factori¬ 
zation  the  ratio  of  the  maximum  and  minimum  elements  of  D  is  examined. 
For  a  2  x  2  block  in  D  the  square  root  of  the  determinant  of  the  block  is 
used.  It  has  been  found  that  if  this  ratio  exceeds  20000  when  using  a  32 
bit  computer  word  length,  then  the  factorization  may  be  unreliable.  The 
matrix  Z  should  be  properly  scaled  and  this  is  usually  done  automatically 
by  the  CHAOS  program. 

In  section  4  the  various  different  wire  junction  treatments 
are  discussed.  Some  of  these  give  rise  to  unsymmetric  matrices.  For  these 
cases  the  Harwell  routine  MA23  {7}  is  used  to  factorize  the  matrix  as 

Z  -  L.U,  ....(34) 

where  L  is  lower  triangular  and  U  is  upper  triangular.  The  reliability  of 
this  factorization  is  checked  by  examining  the  ratio  of  maximum  to 
minimum  elements  of  the  diagonal  of  U.  If  this  exceeds  20000  using  a  32 
bit  computer  word  length,  then  the  factorization  may  be  unreliable.  Again 
the  matrix  Z  should  be  properly  scaled. 

In  those  cases  where  the  ratio  exceeds  the  value  20000,  the 
calculation  is  immediately  terminated  and  nothing  can  be  done  to  get 
answers,  except  by  modifying  the  computer  program.  This  latter  is  not 
advised. 
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Calculation  of  electric  fields 


Once  the  current  amplitudes  1^  of  the  expansion  (28)  are  known, 
we  can  find  the  charge  distribution  from  equation  (19)  and  the  scattered 
electric  fields  at  any  point  in  space  from  equations  (22),  (23)  and  (24). 
The  calculation  of  fields  far  from  the  scatterer  leads  to  some  simplifi¬ 
cations  in  these  expressions.  From  figure  2  we  see  that  the  distance  R, 
which  is  the  distance  from  a  point  on  the  wire  scatterer  to  the 
measurement  point  (receiver),  is  just 


R2  =  r2  *■  r2  -  2r  r  cos  £  , 
r  n  r  n  ^n 


(35) 


where  r  and  r  are  the  distances  from  the  origin  to  the  point  n  of  the 
— n  — r 

scatterer  and  to  the  receiver  respectively  while  £n  is  the  angle 
subtended  by  these  two  points  at  the  origin.  In  radar  cross-section 
studies,  for  which  CHAOS  is  designed,  the  distances  R  and  rr  will  be  large 
so  that  we  may  expand  R, keeping  terms  to  1/r^  only,  to  get  from  equations 
(23)  and  (28) 


•ikr 


s-  .  e 
A  (r  )  *  y 


I  I.  f  T.U)elkr  COS  *  dl,  ....(36) 
i*i  1  -----  -1 
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4>S(r  )  *  o(~) .  ....(37) 
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Thus,  the  components  of  the  far  electric  field  in  a  spherical  coordinate 
system  become 


K  '  »(±t) .  > 


r  '  2 
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....(38 


The  transverse  component  of  the  electric  field  at  large  distances  is 
-ikr 


4itr 


I  ?i.  !  VT.c'-ikr  co* 5 


e\u  -  -  io.y  l  I;  J  u,.T.  ww"  ”  dl,  ....(39 


r  i=*i  wires 


where  ur  is  some  unit  vector  transverse  to  rf,  the  direction  of  the 
receiver. 
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The  program  could  be  modified  to  give  near  field*  a*  well  a* 
far  but  this  has  not  been  done.  Should  it  ever  be  attempted  in  future, 
the  fields  close  to  the  structure  will  exhibit  singularities  and  may  be 
unreliable  since  the  charge  distribution  on  the  wire  structure  goes 
through  a  discontinuous  change  at  the  start  and  end  of  each  triangle 
function.  It  would  be  necessary  for  very  careful  thought  to  be  given  to 
the  problem  of  avoiding  wrong  results  due  to  these  discontinuities.  It 
may  only  be  possible  to  obtain  fields  on  a  grid  of  carefully  chosen 
points  close  to  the  structure  with  interpolation  being  used  for  points 
not  on  the  grid. 

3.6  Radar  cross-sections 


The  radar  cross-section  g(6,$)  is  defined  as  the  area  for  which 
the  incident  wave  contains  sufficient  power  to  produce,  by  omni¬ 
directional  radiation,  the  same  back  scattered  power  density.  So  we  have 


W  I Es.u  |a 
g(e,<j>)  -  — -  "*  -r' 


..(40) 


P.  » 
in 


where  the  power  of  the  incident  wave  P.  is 
r  in 


P. 

in 


IV* 


.(41) 


*  —  for  a  unit  plane  wave, 


where  n  ■  ....(42) 

CHAOS  will  accept  any  polarisation,  linear,  elliptic  or  circular,  for  the 
incident  plane  wave  and  will  calculate  the  radar  cross-section  for  a 
receiver  having  any  linear,  elliptic  or  circular  polarisation.  The 
receiver  need  not  be  in  the  same  placet  as  the  transmitter  so  that  CHAOS 
will  calculate  bistatic  as  well  as  monostatic  radar  cross-sections. 


3.7  Numerical  evaluation  of  impedance  matrix  elements 

The  numerical  evaluation  of  the  impedance  teatrix  elements 

given  by  equation  (32)  needs  some  explanation  since  the  integrand 

contains  a  singularity  and  involves  a  double  integration.  Let  us  write 

Z. .  as 
Ji 
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wires  wires 


Z*?.  -  f  dZZ  .T.(t).T.U). 
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.... (46) 


The  three  terms  are  due  to  the  magnetic  vector  potential,  the  electric 
scalar  potential  and  finite  conductivity  respectively  which  accounts  for 
the  superfixes  of  M,  E  and  C.  As  each  triangle  function  T.  is  non-zero 
only  over  four  segments  (see  figure  1(a)),  the  integrations  reduce  to 
integrals  over  at  most  eight  different  segments.  These  are  replaced  by 
sums  over  these  segments  by  representing  the  basis  function  and  its 
derivative  by  a  sum  of  four  step  functions,  as  illustrated  in  figure  3. 
The  weight  function  and  its  derivative  are  replaced  by  a  sum  of  four 
impulses,  as  shown  in  figure  4.  This  enables  the  double  integrals  in 
equations  (44)  and  (45)  to  be  replaced  by  a  sum  of  sixteen  terms  as 
follows  s- 


z”  *  wop  l  l  t.  t.  D*Vn, 
J1  m»i  n-i  J“  in  J°  JB 

ZE.  *  "T““  l  l  t;  t:  \|Kn , 

J1  wt  «£,  n£,  J® 


....(47) 


....(48) 


where  D.  is  the  scalar  product  between  the  vectors  in  the  directions 
jm 

r-  -  r.  and  r.  -  r.  and  where  r.  is  the  vector  from  the 
±-j,m+i  —  j,m  — i,n+i  —in  — jm 

origin  to  the  start  of  segment  m  under  triangle  function  j.  This  notation 
is  illustrated  in  figure  5.  Finally,  the  i|>  function  is  defined  as 


i  in 

*jm  "  4*7 


"2 i-i +n 


_  r^2 i-J  +n 

*'2i-2+n>  Si-.+n 


in  which  R  is  the  distance  between  the  field  point  adjacent  to  the 
midpoint  of  segment  m  under  triangle  function  j  and  the  source  point  on 
the  nc  segment  under  triangle  function  i.  The  index  m  causes  the 
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•lunation  to  extend  over  the  four  segments  forming  the  baae  of  triangle 
function  j  while  the  index  n  causes  the  summation  to  extend  over  the  four 
segments  forming  the  base  of  triangle  function  i.  The  formulae  used  to 
evaluate  the  integral  in  equation  (49)  have  been  derived  by 
Harrington  [8]  and  are  tabulated  in  appendix  A. 

The  alert  reader  will  have  noticed  the  use  of  the  rather 
imprecise  term  "adjacent"  in  the  previous  paragraph.  The  integral  in 
equation  (49)  contains  a  singularity  when  the  field  point  coincides  with 
a  point  on  the  line  source.  In  order  to  remove  this  singularity,  it  has 
generally  been  decided  to  consider  the  line  sources  as  situated  on  the 
axes  of  the  wires  while  the  field  points  are  situated  on  the  surfaces. 
The  position  of  the  field  point  on  the  wire  surface  is  illustrated  in 
figure  6  and  is  determined  as  follows:  a  line  is  drawn  from  the  source 
point  to  the  axis  of  the  wire  on  whose  surface  the  field  point  is  to  be 
situated.  At  right-angles  to  both  this  line  and  this  axis,  another  line 
is  drawn.  This  line  cuts  the  wire  surface  in  two  points.  One  of  these 
points  is  taken  as  the  field  point;  while  there  is  still  ambiguity  in 
which  point  is  the  field  point,  there  is  no  ambiguity  in  the  distance  R. 
I"  is  the  same  from  both  points.  In  the  degenerate  case  of  a  field  point 
on  the  same  segment  as  the  source,  any  point  on  the  circumference  will  do 
and  each  point  gives  the  same  distance  R. 

For  two  reasons  the  impedance  matrix  calculated  in  the  manner 
described  here  will  no  longer  be  exactly  symmetric.  The  first  is  that  the 
mathematical  representation  of  the  weighting  function  (shown  in 
figure  4)  is  different  from  that  used  for  the  basis  function  (shown  in 
figure  3).  The  second  is  the  method  of  avoiding  the  singularity  in 
equation  (49).  Harrington  [2]  has  found  that  little  difference  is  made  to 
the  radar  cross-sections  by  averaging  the  off-diagonal  elements  of  the 
impedance  matrix.  This  has  also  been  confirmed  by  our  own  calculations 
and  this  averaging  is  usually  performed  by  CHAOS  although  the  unsymmetric 
matrix  can  be  used  if  desired  by  specifying  the  appropriate  option.  The 
advantages  of  symmetry  are  substantially  reduced  computer  storage  needs 
and  the  exact  satisfaction  of  the  principle  of  reciprocity. 
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3.8  Numerical  evaluation  of  the  voltage  vector 

The  numerical  evaluation  of  the  voltage  vector  given  by 
equation  (31)  ia  straightforward.  In  the  CHAOS  program  the  source  is 
assumed  to  be  a  plane  wave  of  the  form 

E1(r)  -  E  e'1-*-.  _ (50) 


The  form  of  ^  determines  the  polarisation  of  the  wave  as  described  in 
appendix  B.  The  weight  function  is  represented  by  a  set  of  impulse 
functions  (as  in  figure  4)  and  so  we  get 
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where  A.  is  the  vector  representing  the  length  and  direction  of  the  m1*1 
— jm 

segment  of  triangle  function  j. 


4.  JUNCTIONS  OF  WIRES 


There  are  two  conditions  that  must  be  satisfied  at  a  junction. 
The  first  is  Kirchoff's  law  which  states  that  the  sum  of  the  currents 
entering  a  junction  equals  the  sum  of  the  currents  leaving  the  junction. 
The  second  states  that  the  scalar  potential  is  continuous  across  the  wire 
surfaces  at  the  junction.  Kirchoff's  law  is  relatively  easy  to  satisfy 
but  the  condition  on  scalar  potential  has  not  been  satisfied  in  any 
computer  program  known  to  the  author.  This  is  because  the  unknown 
variables  in  method  of  moments  solutions  are  the  currents  and  their 
derivatives,  the  charges.  Thus,  we  really  need  a  condition  on  the  charge 
not  the  scalar  potential.  It  is  possible  to  write  down  an  analytical 
expression  for  the  charge  but  it  is  not  very  suitable  for  a  computer 
program.  So  various  authors  have  put  forward  a  variety  of  approxi¬ 
mations. 


Two  in  particular  have  been  picked  out  and  programmed  into 
CHAOS  but  neither  is  entirely  satisfactory  so  an  approximation  based  on  a 
combination  of  these  two  methods  has  also  been  programmed.  Three  other 
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ideas  have  also  been  tried  out  but  are  only  of  limited  use.  In  none  of 
these  methods  is  it  necessary  to  place  further  restrictions  on  the  length 
of  segment  or  radius  of  the  wire  abutting  the  junction  than  those  which 
must  be  satisfied  at  any  point  of  a  wire,  as  described  in  appendix  A. 

4.1  Chao  and  Strait's  method 


The  first  approximation  is  used  by  both  Chao  and  Strait  [1]  and 
Richmond  [4],  Suppose  N  wires  meet  at  some  junction,  then  Kirchoff's  law 
is  satisfied  exactly  by  overlapping  wire  i  by  two  segments  on  to 
wire  i  -  1,  except  that  wire  1  is  not  overlapped  on  to  wire  N.  This  is 
shown  in  figure  7(a)  for  three  wires.  This  allows  extra  triangle 
functions  on  wires  2  and  3  which  straddle  the  junction  (these  are 
labelled  I  and  IT  in  figure  7(a)).  The  total  current  l(Jl)  on  a  segment 
adjacent  to  a  junction  is  the  sum  of  the  current  on  that  segment  and  the 
extra  overlap  segment  if  it  exists.  The  line  density  of  charge  is 
obtained  from  equation  (19)  which  is 


a(£) 


1  dI<*> 

id  di 


As  the  wire  connection  data  are  read  in,  the  program  records  in 
an  array  the  wire  numbers  of  all  wires  that  enter  each  junction  in  the 
order  in  which  they  are  encountered  in  the  input  data.  This  is  slightly 
modified  for  problems  with  rotational  symmetry  (see  section  7). 

An  aspect  of  this  overlapping  wire  method  that  does  not  appeal 
to  a  mathematician  is  the  asymmetry  and  arbitrariness  introduced  by 
having  to  choose  one  wire  at  each  junction  to  have  no  overlap  and  by  the 
choice  of  which  order  the  remaining  wires  overlap  each  other.  Two 
different  people  setting  up  the  same  problem  will  not  necessarily  number 
the  wires  and  junctions  in  the  same  way  so  that  a  different  set  of 
overlaps  will  result  and  a  different  set  of  triangle  functions  will  be 
set  up.  But  the  sum  of  the  appropriate  two  triangle  functions  on  a 
segment  adjacent  to  a  junction  will  Btill  give  the  same  total  current  - 
at  least  to  within  the  accuracy  of  the  calculation.  This  objection  has  a 
practical  effect  when  dealing  with  a  structure  having  left-right  symmetry 
or  N-fold  rotational  symmetry  which  is  discussed  further  in  section  7. 
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The  only  satisfactory  way  of  removing  this  ambiguity  is  to  do 
away  with  overlaps  completely.  However,  this  means  an  extra  constraint 
(representing  Kirchoff's  law)  must  be  introduced  and  this  leads  to  an 
extra  row  in  the  impedance  matrix  whose  elements  may  be  of  very  different 
magnitudes  since  they  do  not  represent  impedances.  Also  for  a  triangle 
function  basis,  one  must  introduce  right-angled  triangle  functions  at 
the  ends  of  each  wire  which  complicates  the  programming.  Thus,  in  CHAOS 
we  generally  use  overlapping  triangle  functions  and  live  with  the 
arbitrariness.  These  overlaps  are  worked  out  automatically  by  the 
geometry  package  in  the  CHAOS  program  and  the  program  user  need  not  worry 
over  this. 

4.2  Sayre's  method 

The  second  method  has  been  suggested  by  both  Curtis  [9]  and 
Sayre  [10]  while  an  improved  scheme  has  been  put  forward  by  King  [11]. 
These  methods  all  involve  charge  redistribution.  One  starts  from 
equation  (19)  and  at  a  junction  of  N  wires  integrates  this  over  a  small 
volume  surrounding  the  junction.  This  volume  is  defined  by  the  ends  of 
the  N  segments  abutting  the  junction.  Sayre  makes  the  assumption  that  the 
line  density  of  charge  is  constant  over  these  N  segments  and  arrives  at 
the  expression 
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where  £  .  is  the  position  of  the  other  end  of  the  segment  abutting  the 
junction  on  wire  uk  and  A  .  is  the  length  of  this  segment.  Now  because  of 
the  overlaps  used  to  satisfy  Kirchoff's  law,  there  would  appear  to  be 
2N  -  1  segments  abutting  a  junction  of  N  wires.  For  the  purpose  of 
representing  this  average  charge,  the  N  -  1  segments  due  to  overlaps  are 
ignored  and  no  charge  is  spread  over  them.  The  charge  is  then  spread 
uniformly  over  the  N  segments  adjacent  to  the  junction  which  were  on  the 
original  wires. 
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Curtis' 8  scheme  made  the  assumption  that  the  surface  density 
of  charge  was  constant,  while  King,  working  with  a  tapered  wire,  produces 
two  formulae,  one  for  short  wires  the  other  for  long  ones.  In  the  limit 
of  wires  of  equal  radius  King's  formula  reduces  to  Sayre's  showing  that 
there  is  some  justification  for  assuming  that  the  line  density  of  charge 
is  constant.  However,  King's  scheme  has  not  been  programmed  because  it  is 
known  that  Sayre's  treatment  does  not  always  give  acceptable  results  for 
junctions  of  wires  of  equal  radius. 

4.3  50/50  method 


Since  neither  Chao  and  Strait's  nor  Sayre's  methods  always 
give  acceptable  results,  a  new  scheme  was  tried.  It  has  been  noticed  that 
when  Chao  and  Strait's  method  disagrees  with  experiment,  Sayre's  method 
usually  gives  acceptable  results.  The  converse  also  appears  to  hold. 
Thus,  it  was  decided  to  calculate  the  charge  by  both  methods  and  combine 
them  with  equal  weight  to  give  the  formula 

\  oU)  +  i  a.  ....(53) 

Thus,  in  calculating  the  charge,  one  is  taking  account  of  both  the 
individuality  of  the  wire  and  of  the  aggregate  of  wires  meeting  at  the 
junction.  No  attempts  have  been  made  to  optimize  the  expression  (53)  by 
choosing  different  weightings  (eg,  30/70,  80/20,  etc). 

4.4  Other  methods 


It  has  been  noticed  that  the  Chao  and  Strait  method  is 
generally  better  for  wire  meshes  representing  solid  surfaces  while  the 
50/50  method  is  generally  better  for  junctions  of  actual  wires.  For  these 
reasons  an  option  is  provided  in  CHAOS  which  allows  certain  junctions  to 
use  Chao  and  Strait's  methods  and  others  to  use  the  50/50  method. 

In  another  scheme  the  currents  on  all  wires  entering  a  junction 
are  forced  to  be  zero.  This  means  that  the  overlapping  segments  are  not 
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required.  Sayre's  method  for  charge  redistribution  is  used. 
Surprisingly,  this  scheme  gives  good  results  for  some  geometries; 
however,  its  use  is  not  recommended. 

Finally,  a  more  promising  method  without  overlaps  is 
programmed.  In  this  scheme,  each  triangle  function  abutting  a  junction 
of  two  or  more  wires  was  replaced  by  a  step  function  on  the  two  segments 
adjacent  to  the  junction,  as  shown  in  figure  7(b),  while  Sayre's  charge 
redistribution  scheme  is  used.  The  saving  in  matrix  size  gained  by  the 
simpler  junction  treatment  was  offset  by  the  need  for  more  segments  along 
each  wire.  The  method  required  6  and  preferably  8  segments  between  any 
two  junctions  whereas  the  50/50  method  only  needed  2  on  the  same  problem. 
Also  Kirchoff's  law  is  not  satisfied  exactly  so,  for  these  reasons,  the 
method  has  not  been  pursued  but  it  has  been  left  in  the  CHAOS  program. 


4.5 


Effects  on  impedance  matrix  elements 


As  the  Chao  and  Strait  junction  method  does  not  involve  charge 
re-distribution,  it  does  not  affect  the  way  the  impedance  matrix  elements 
are  calculated,  but  all  the  other  methods  do.  Of  the  three  terms  that  the 
impedance  matrix  elements  have  been  split  into,  it  is  only  that  due  to 
the  electric  scalar  potential  that  is  affected.  The  factor  (d/dd' )T^(d' ) 
in  equation  (45)  and  t!  in  equation  (48)  need  to  be  replaced.  In  order 
to  see  the  form  these  take  let  us  suppose  that  at  some  junction  of  N 
wires,  triangle  functions  Jj  to  JM  are  directed  out  of  the  junction  (ie, 

Ni 

leave  the  junction),  Kj  to  K^_1  straddle  the  junction  (ie,  are  overlaps) 

and  L.  to  L  are  directed  into  the  junction  (ie,  enter  the  junction).  As 
n2 

N  wires  meet  at  the  junction,  Ni  +  N2  -  N.  Now  the  current  leaving  the 
volume  surrounding  the  junction  (as  defined  in  section  4.2)  is  given  by 
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From  the  form  of  equation  (54)  it  is  clear  that  the  modification  will 
affect  2N  -  1  columns  of  the  impedance  matrix  since  for  Sayre's  method 
equation  (45)  becomes 
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where  I  is  given  by  equation  (54).  The  a(  i  ')  and  b(  2  ')  are  defined  by:- 


(  1  normally, 
a(2 ' )  =  ( 

(0  if  2'  is  on  a  segment  adjacent  to  a  junction, 

(  0  normally, 

b(r)  =  ( 

(1  if  2'  is  on  a  wire  segment  (not  an  overlap  segment) 
adjacent  to  a  junction. 


This  formula,  equation  (55),  is  not  quite  general  since  a  triangle 
function  can  leave,  straddle  and  enter  three  junctions  (see  figure  8) 
although  only  two  junctions  will  contribute  to  the  average  charge  since 
we  have  decided  not  to  spread  any  charge  on  an  overlap  segment  (section 
3.2).  So  a  further  term  in  the  curly  brackets  is  required  similar  to  the 
b(2')  term  but  referring  to  this  second  junction. 

The  equation  corresponding  to  equation  (48)  is 


<» 


Z?.I.  =  lit-  (a.(n)t!  I.  +  b.(n)Q}\(i.  , 

ji  i  iU)e  im  i  in  i  i  jm 

J  m=i  n=i  J  J 


.  (56) 


where 


(2,  .  -  2  .  )  r  N!  N2  n 

q  =  -ii-iiy-  r  r  t:  I.  +  i  t:  1.1 . 

N  L-  T  i1  J  •  T  I1*  1J 

,  V  A  LJ=J1  3  3  J=L1  J 


2  l  L 

P=1 up 


....(57) 
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The  terms  involving  the  K  type  triangle  functions  (ie,  those  that 
straddle  the  junction)  are  identically  zero  because  the  program,  as  it  is 
at  present  set  up,  makes  all  segments  between  any  two  junctions  equal  in 
length.  Should  this  be  changed  in  future,  then  an  extra  term  will  need  to 
be  added  into  equation  (57).  From  figure  9  we  see  that:- 

a.(n)  =  1  normally, 

=0  if  triangle  function  i  leaves  a  junction  J  and  n  =  1 
(cases  (i)  and  ( i i ) ) , 

=0  if  triangle  function  i  straddles  a  junction  K  and 
n  =  2  or  3  (cases  (iii)  and  (iv)), 

=0  if  triangle  function  i  enters  a  junction  L  and  n  *  4 
(cases  (v)  and  (vi)), 

and  :- 

b^(n)  =  0  normally, 

=1  if  triangle  function  i  leaves  a  junction  J  and  if 
wire  on  which  i  is  situated  has  no  overlap  at  J 
and  if  n  -  1  (case  ( i ) ) , 

=1  if  triangle  function  i  straddles  a  junction  K 

and  either  n  =  2  if  wire  enters  junction  (case  (iii)) 
or  n  =  3  if  wire  leaves  junction  (case  (iv)), 

=1  if  triangle  function  i  enters  a  junction  L  and  if 
wire  on  which  i  is  situated  has  no  overlap  at  L 
and  if  n  =  4  (case  (vi)). 

In  retrospect,  the  situations  for  which  b^(n)  are  non-zero  could  be  more 
logically  chosen.  One  might  argue  that  in  case  (ii)  b^  (1)  =  1,  while  in 
case  (iv)  b^  (3)  =  0,  and  similarly  in  case  (v)  b^  (4)  =  1,  with  case 

(iii)  b^  (2)  =  0.  Alternatively,  one  might  argue  that  in  cases  (ii)  and 

(iv)  b^  (1)  =  p  and  b^  (3)  =  1  -  p  respectively  with  a  similar  split  for 
cases  (v)  and  (iii).  One  would  perhaps  like  to  choose  this  weight  p  as 
the  ratio  of  current  magnitude  due  to  one  triangle  function  to  the  sum  of 
the  magnitudes  of  the  two  triangle  functions  which  are  based  on  the 
segment  adjacent  to  the  junction.  However,  this  produces  a  non-linear 
problem  and  is  to  be  avoided.  For  better  or  for  worse,  we  have  chosen 
this  weight  p  to  be  zero  and  no  study  has  been  made  of  other 
alternatives;  perhaps  one  day  this  should  be  done. 

As  with  equation  (55),  expression  (56)  is  not  quite  general 
since  two  of  the  situations  listed  above  could  occur  simultaneously.  This 
could  arise  since  we  know  from  figure  8  that  a  triangle  function  can 
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leave  a  junction  J,  straddle  a  junction  K  and  enter  a  junction  L.  As 
depicted  there  would  be  contributions  to  triangle  functions  i  from  both 
junctions  J  and  K  but  not  L.  The  converse  could  also  happen  in  which 
contributions  from  junctions  L  and  K  occurred  but  not  J.  The  third 
possibility  arises  when  a  triangle  function  leaves  a  junction  J,  does  not 
straddle  any  junction  but  enters  another  junction  L.  Thus,  another  term 


of  the  type  b^(n)Q  referring  to  this  second  junction  needs  to  be  included 


in  the  curly  brackets  of  equation  (56). 


However,  the  50/50  method  produces  a  further  complication; 
expressions  (56)  and  (57)  remain  the  same  but  the  coefficients  a^(n)  and 
b^(n)  are  modified.  Where  the  parameter  a^(n)  *  0  in  the  Sayre 

formulation  it  is  set  to  0.5  in  the  50/50  method.  Similarly,  where  the 


parameter  b.(n)  *  1  in  the  Sayre  formulation  it  is  set  to  0.5  in  the 


50/50  method. 


Finally,  it  is  clear  from  equation  (56)  that  the  resulting 


matrix  Z.^  will  be  unsymmetric  since  the  original  symmetric  matrix  has 
been  replaced  by  a  matrix  whose  columns  are  linear  combinations  of  the 


original  columns  without  any  corresponding  replacement  of  rows.  The 
Harwell  matrix  factorization  package  MA23  is  used  for  this  situation. 


More  recently  tests  have  suggested  that  the  result  of  averaging  the  off 


,E 


diagonal  terms  to  force  the  matrix  Z^  to  be  symmetric  is  in  no  worse 


agreement  with  experiment  and  in  some  ways  is  in  rather  better  agreement 
(in  the  sense  that  reciprocity  is  satisfied).  Thus,  the  option  exists  to 
force  the  matrix  to  be  symmetric  and,  in  view  of  the  saving  in  computer 
storage  that  follows  from  this,  it  is  now  the  recommended  option. 


WIRE  MESH  MODELLING  OF  SOLID  SURFACES 


On  the  face  of  it  a  computer  program  based  on  the  thin  wire 
approximation  is  inappropriate  for  calculations  involving  surfaces. 
However,  it  is  believed  that  at  present  no  reliable  general  purpose 
surface  patch  computer  code  has  been  written  which  can  cope  with  finite 
conductivity  or  which  results  in  a  symmetric  impedance  matrix.  The  reason 
for  this  is  essentially  that  the  calculation  of  the  elements  of  the 


impedance  matrix  would  be  very  expensive  and  perhaps  difficult  as  the 
expression  analogous  to  equation  (32)  involves  a  four-fold  rather  than  a 
two-fold  integral.  Intuitively,  one  feels  that  a  smaller,  perhaps  much 
smaller,  matrix  would  result  so  that  for  problems  involving  many  aspect 
angles  and  polarisation  combinations,  a  surface  patch  approach  might  be 
viable. 


So  how  does  one  use  CHAOS  to  model  a  surface  of  a  conductor 
which  is  not  necessarily  a  perfect  conductor?  The  trick  is  to  recall  that 
a  wire  mesh  in  the  same  shape  as  a  solid  surface  will  give  very  similar 
radar  cross-sections  provided  the  mesh  size  is  small  compared  with  the 
wavelength.  Thus,  one  replaces  the  solid  surface  by  a  wire  mesh  whose 
holes  are  roughly  0.02  to  0.03  square  wavelengths  in  area.  This  size  has 
been  determined  by  experience.  The  main  worry  with  this  approach  is  that 
one  is  forcing  the  currents  to  flow  along  the  wires  of  the  mesh  rather 
than  as  they  would  in  real  life.  For  this  reason  a  triangular  mesh  rather 
than  a  rectangular  one  is  preferred.  Also  one  would  tend  to  use  a  finer 
mesh  on  surfaces  of  high  curvature  than  on  flat  surfaces,  and  at  edges  of 
surfaces  where  larger  currents  might  be  expected  to  flow.  Should  any  part 
of  the  surface  suggest  a  preferred  direction  for  the  current  to  flow, 
then  the  mesh  sides  should  be  aligned  in  this  direction. 

Next,  there  is  the  question  of  what  radius  to  give  the  wires  in 
the  mesh.  Experience  has  again  shown  that  this  should  be  chosen  so  that 
the  surface  area  of  the  wire  mesh  is  no  less  than  the  area  of  the  solid 
surface.  In  this  way  the  position  of  the  lobes  of  the  radar  cross- 
section  pattern  will  be  correctly  predicted.  However,  only  the  amplitude 
of  the  major  lobes  will  be  reliably  calculated.  If  the  radius  has  not 
been  correctly  chosen,  one  may  expect  errors  in  amplitude  of  a  few 
decibels  for  lobes  10  db  below  the  maximum  while  for  smaller  lobes,  say 
20  to  30  db  below  the  maximum,  errors  of  5  to  10  db  have  been  seen.  In 
many  applications  this  sort  of  accuracy  is  quite  acceptable.  So  no  exact 
prescription  can  be  given  for  the  radius  of  the  wires  in  the  mesh  and  up 
to  five  times  that  suggested  above  may  be  necessary  to  get  good  agreement 
with  experiment  on  the  minor  lobe  amplitudes.  This  sort  of  variation  in 
wire  mesh  radius  has  very  little  effect  on  the  amplitudes  of  the  major 
lobes,  perhaps  of  the  order  of  a  decibel. 
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The  wires  should  be  given  the  seme  conductivity  ss  the  solid 
surface.  In  the  case  where  the  surface  is  composed  of  strips  of  material 
having  different  conductivities  along  and  across  the  strips,  then  the 
mesh  should  be  rectangular  in  shape  aligned  along  and  across  the  strips, 
the  wires  in  one  direction  having  one  conductivity,  the  ones  in  the  other 
the  other  conductivity. 

Finally,  we  have  found  that  the  Chao  and  Strait  junction 
treatment  for  wires  in  meshes  representing  solid  surfaces  tends  to  give 
much  better  predictions  for  minor  lobes  than  does  Sayre's  treatment. 

6.  h  TREATMENT  OF  MATERIALS  WHICH  ARE  NOT  PERFECTLY  CONDUCTING 

6.1  Types  of  current 

At  this  point  one  must  consider  the  various  types  of  complex 
current  density.  Inside  a  conductor,  the  current  is  almost  entirely  due 
to  the  motion  of  free  electrons  and  is  called  the  conduction  current  and 
is  given  by 

J  *  oE,  ....(58) 

where  o  is  the  conductivity. 

In  free  space  we  have  no  motion  of  charges  at  all  and  we  have 
only  a  "free  space  displacement  current"  given  by 

iwe  E,  ....(59) 

o— 

where  e  is  the  dielectric  constant  of  free  space  (»  (l/36ir)  x  10-9 
o 

F/m) . 


In  matter,  in  addition  to  the  conduction  current  and  the  free 
space  displacement  current,  there  is  a  current  due  to  the  motion  of  the 
bound  charges.  This  is  called  the  polarisation  current  and  is  given  by 


iw(e  -  e  )E. 

O  “ 


(60) 
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For  completeness  one  should  note  that  the  current  itsE  is  called  the 
displacement  current. 

Both  the  conduction  and  polarisation  currents  give  rise  to 
scattered  electric  fields  due  to  the  motion  of  the  charges.  So  for  a  non¬ 
magnetic  material  we  would  like  to  calculate  both  these  currents.  Let  us 
write 


J8  *  J  ♦  ito(e  -  cq)E»  ....(61) 

g 

then  J  gives  rise  to  the  total  scattered  field  for  a  nonmagnetic 
material . 


6.2  Impedance  of  wire 

By  substituting  for  from  equation  (58)  into  the  third  of 
equations  (6)  and  also  using  the  first  we  have 

V.V.E  -  T2E  *  0,  - (62) 

where  T2  *  -  itop(o  +  iwe  )  •  Clearly,  we  also  have 

-  T2J  *0  •  •  • .(63) 

and  so  V„V„j8  -  T2J8  =  0.  ....(64) 

Now  let  us  specialize  to  a  cylindrical  wire  of  infinite  length  lying 
along  the  Z  axis.  Let  us  assume  the  current  only  flows  in  the  Z  direction 
and  is  a  function  of  p,  the  radial  distance  only.  Let  us  use  cylindrical 
polar  coordinates  (p,$,Z)  based  on  the  Z  axis.  Then 

J8  *  (0,0, J(p))  - (65) 

and  equation  (64)  reduces  to 

?  ip (p  $  *  T’J  ’  »•  ••••«*> 


29 


This  is  Bessel's  equation.  We  require  that  solution  which  is  finite  at 

the  centre  and  which  gives  the  current  density  i  ,  say,  at  the  surface  of 

o 

the  wire.  So 

J  <Tp) 
o 

where  a  is  the  radius  of  the  wire  and  Jq  is  the  cylindrical  Bessel 
function  of  order  zero  which  is  regular  at  the  origin. 


The  total  current,  I,  flowing  in  the  wire  is  given  by 


r2lT  ,a 

I  *  J #  d$  J#  pdpj(p) , 


•  2ira 
Xo  T 


Jj(Ta) 

J  (fa7’ 
o 


....(68) 


where  Ji  is  the  cylindrical  Bessel  function  of  order  unity  which  is 
regular  at  the  origin.  The  electric  field  at  the  wire  surface  is  given  by 


h> 


(o  + 


iw(e  -  €. 


7T 


V» 


(69) 


where  the  wire  impedance  per  unit  length  Z  is  given  by 


u 


J0(T«) 


0] 


2Trafo  +  iw(e  -  e  ))  J,  (Ta)’ 

O  * 


....(70) 


It  has  been  found  that  this  formula  is  inadequate  for  some  materials,  in 
particular  those  which  are  not  solid  but  consist  of  strands  of  fibre 
woven  together.  For  these  materials  the  amended  formula 


t  J„(Ta) 

_ £ _ O 

2ira(a  +  iw(e  -  e  ))  •  J.<Ta> 
o  * 


(l  + 


(71) 


is  used  where  v  is  the  "frequency  fit"  parameter  in  cycles  per  second  and 
is  chosen  to  give  optimum  agreement  with  experimental  impedance  per  unit 
length  data. 
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The  use  of  this  formula  has  been  found  to  give  good  radar 

cross-sections  for  "wires"  ranging  from  copper  (a  ■  5.8  x  107  mho/m)  to 

salt  water  (o  *  20.37  mho/m,  e  *  50e  )  and  polystyrene  (e  *  2.54e  ,  loss 

o  o 

tangent  =  4.05  x  10  "). 

No  computer  routines  have  been  found  that  will  accurately 
calculate  the  Bessel  functions  JQ  and  Ji  for  an  arbitrary  complex 
argument,  £;  however  Che  ratio  Jj/Jo  can  be  found  satisfactorily  for  all 
arguments  from  the  continued  fraction  (12) 


J 


J  (c) 

o', 

(0 
q-l  s 


C 

2q  -  C2 

2(q+l)  -  C2 

2(q+2>  -  C2 

2(q+»)  -  £2. 


(72) 


For  large 


Id  we  have 


_  cos  (c  -  fff)  + 

j  Co  cob  (r,  -  trri 

o 


o(z-1), 


a  i  for  large  Im(^)  >  0, 

=  -i  for  large  Im(^)  <  0. 


see* 


(73) 


For  small  |^j  we  have 

•>,'<>  i 


•  •  •  • 


(74) 


Some  other  quantities  are  perhaps  worth  recording.  The  skin 
depth,  6,  is  given  by 


6  = 


e  •  a  a 


(75) 


The  loss  tangent  is  the  ratio  of  imaginary  to  real  components 
of  the  dielectric  constant  and  the  loss  angle  n  is  given  by 


tan  n 


Im(e) 


•  t  t  • 


(76) 
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The  resistance  Ry  and  inductance  per  unit  length  are  related  t°  Zw  by 

\  +  i“Lw  "  Zu)*  ....(77) 

The  dc  resistance  and  inductance  per  unit  length  are  given  by 


V  "  2 

....(78) 

iroa 

ldc  -  fcr  <“  50  nH/n)* 

....(79) 

At  high  frequency  and  provided  a  »  U)e,  we  have  approximately 


6.3 


Form  of  matrix  of  contributions  due  to  finite  conductivity 


(80) 


The  matrix  of  contributions  due  to  finite  conductivity  is 
given  by  substituting  equation  (71)  into  equation  (46).  Because  the 
current  basis  functions  are  a  set  of  overlapping  triangles,  this  matrix 
is  sparse  with  most  contributions  on  or  immediately  adjacent  to  the  main 
diagonal.  In  the  special  case  of  a  dipole  the  matrix  is  tridiagonal,  but 
for  problems  involving  three  or  more  wires  joined  at  a  vertex  this  ceases 
to  be  so  and  contributions  outside  the  diagonal  band  occur.  Since  the 
programming  required  to  calculate  correctly  these  contributions  is 
somewhat  involved,  it  would  have  been  easier  to  treat  finite  conductivity 
as  a  series  of  lumped  loads  situated  at  the  centres  of  each  triangle 
function  which  leads  to  a  nearly  diagonal  matrix.  Although  this  has  been 
done  in  other  computer  programs,  this  slightly  less  realistic  approach 
has  not  been  adopted  in  CHAOS. 


6.4  Lumped  loads 


Some  problems  involve  the  attachment  of  loads  to  the 
scattering  body  and  so  an  option  is  provided  to  allow  these  to  be 
included.  These  might  take  the  form 


Z 


L 


R  +  iuL  - 


1 

ituC’ 


*  •  •  • 


(81) 
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where  R,  L  and  C  are  respectively  resistance,  inductance  and  capacitance. 
In  CHAOS  such  loads  may  be  attached  either  at  the  beginning  of  a  wire  or 
at  the  end  of  an  even  numbered  segment.  This  means  that  the  load  is 
attached  in  series  between  the  two  segments.  A  load  attached  at  a  free 
end  of  a  wire  will  have  no  effect  since  the  current  is  assumed  to  be  zero 
at  free  ends. 

The  effect  on  the  impedance  matrix  is  to  alter  the  diagonal 
term  corresponding  to  the  triangle  function  at  whose  centre  the  load  is 
attached.  At  wire  junctions  where  three  or  more  wires  meet,  the  load  may 
be  attached  at  the  centre  of  two  triangle  functions  thus  affecting  two 
diagonal  terms  and  an  off  diagonal  term.  For  instance,  in  figure  7(a),  if 
a  load  were  attached  at  the  end  of  wire  2  which  enters  the  junction,  then 
the  S  and  T  diagonal  elements  of  the  impedance  matrix  will  be  affected  as 
well  as  the  ST  off  diagonal  element.  This  is  because  both  triangle 

functions  I„  and  I„  have  a  common  centre  at  the  end  of  wire  2. 

S  T 

Another  unwelcome  feature  of  lumped  loads  is  that  the 
impedance  matrix  can  become  badly  scaled  as  a  result  of  adding  in  these 
loads.  Experience  shows  that  the  impedance  matrix  should  be  properly 
scaled  before  factorization  but  even  this  may  not  be  sufficient  for 
computers  with  a  short  word  length  (eg,  32  bits)  and  errors  of  up  to  0.5 
db  have  been  seen  in  the  peak  values  of  some  lobes.  Thus,  jobs  involving 
loaded  structures  should  ideally  be  run  on  computers  with  a  longer  word 
length  (eg,  64  bits). 

7.  SYMMETRIES  IN  GEOMETRICAL  STRUCTURE 


In  some  problems  the  geometrical  structure  possesses  special 
symmetries  which  can  be  utilised  to  reduce  substantially  computing  time 
and  storage.  We  consider  first  the  case  of  a  body  with  rotational 
symmetry,  eg,  a  cylinder}  then  a  body  with  left-right  symmetry,  eg,  an 
aircraft;  and  finally  a  body  with  two  planes  of  mirror  symmetry,  eg,  an 
aircraft  sitting  on  a  perfectly  conducting  ground  plane. 
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7.1  Rotation*!  lyttry 


Let  us  consider  a  body  with  n-fold  rotational  ayunetry  about 
some  axis.  In  the  CHAOS  program  this  axis  is  required  to  be  the  Z  axis  but 
the  following  analysis  does  not  depend  on  this.  Let  us  divide  the 
structure  up  into  n  sections,  together  with  a  section  describing  those 
triangle  functions  lying  wholly  on  the  axis  of  rotation.  Then  the 
impedance  matrix  may  be  written  in  block  matrix  form  as:- 


r 

ZB 

*1 

zT 

ZA 

ZT 

ZA 

ZA 

Zi  l 

Zl2 

Zin 

ZA 

Z21 

Z22 

Z2„ 

•  •  •  • 
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Z  , 
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where  the  square  matrix  Zfi  represents  the  interactions  between  triangle 

functions  lying  wholly  along  the  axis  of  rotation,  the  square  matrix  Z^ 

represents  the  interactions  between  section  i  and  section  j  of  the 

structure,  and  the  rectangular  matrix  Z^  represents  interactions  between 

the  axis  and  any  one  of  the  n  sections.  The  matrix  Z  is  assumed  to  be 

T 

symmetric  so  that  Z^j  *  where  T  denotes  transpose.  This  means  that  a 
junction  treatment  which  gives  rise  to  a  synmetric  impedance  matrix  must 
be  used.  We  could  have  done  the  analysis  with  an  unsymmetric  matrix  Z  but 
this  requires  roughly  twice  as  much  storage  and  does  not  guarantee  that 
the  principle  of  reciprocity  is  satisfied.  We  can  immediately  write  down 
some  relations  between  the  blocks  Z.j  due  to  the  rotational  symmetry  of 
the  geometry.  These  are 


Z.  . 
ij 


|  Zn+j-i+i 

<  Vi- 


for  i  >  j, 
for  i  <  j, 


....(83) 
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Then  we  see  that 


S.  .  ,  -  S.  S«,  ,  } 

jk  JU’  j 

} 

Sj-l+»,k  ”  SjkSik’  where  *  *  c°®plex  conjugate,  } 

} 
} 
} 

(n  ifj'i,  j 

} 


S .  -  S.  .  , 

jtn,k  J  »k’ 


|  S.  S*  »  { 
k“i  ^  ^  {0  otherwise. 


....(88) 


We  may  now  expand  the  matrix  equation  ZI^  ■  V  to  obtain 

VB  ’  4  (  I  ij)  ' 

J_1  _ (89) 

n 

z  I  ♦  J  Z. . I .  *  V .  x  m  1 9  2 ,  ••••  n • 

A-B  A  ij-J  -l 

We  now  multiply  the  ith  equation  of  this  latter  set  by  Sfc.  and  sum  over  i 
to  obtain 

{ .KiVu  *  j,  j,skizijlj  ■  .KA  k '  *-  2 . ”• 

11  1  1  J  1  1  1  ....(90) 


Wri 


n 

ite  k  =  l  sjiij 

j=i  J  J 


....(91) 


and 


k  *  X  sj*~j’ 

J=1  J  J 


....(92) 


where  J.  and  W.  are  discrete  Fourier  transforms  of  the  currents  ][.  and 

“JC  ^A»  J 


voltages  Vj.  The  inverse  transformations  give 
l  n 

I.  -  -  Y  S*  J„ 

-J  n  L  3i~i 


and 


V.  X  T  s*  W 
“J  n  iSi  jtft* 


....(93) 

....(94) 
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Substituting  for  Jtj  and  Vj  we  get  , 

n<SkiZA-B  +  n  ^  J^ki^i-fc  =  4  k  =  1,  2,  - n,  ..(95) 


where  Q  =  I  Sy ... 

j  =  i 


....(96) 


By  substituting  for  Z. .  and  doing  some  rearrangement  of  the  summation  we 
find 


QU  '  s?t  I  Vp-  ••••<97> 

p-i  r  r 

So  we  want  to  solve  the  equations 

n 

n^ji^A~ B  +  ^  I  **pk4^4  =  — k  k  ~  ^  ....  n  ....(98) 

....(99) 

For  k  =  2,  3  ....  n  we  solve  the  block  diagonal  system  of  equations 

....(100) 

....(101) 


and 


ZbIb  +  ZA— 1  =  — B* 


44  =  4’ 


n 


where  U  =  I  S*  Z  , 

k  pk  p 


while  for  k  =  1  we  solve  the  augmented  equation 

iz  ZT 

n  B  A 

n 

ZA  l  Zt 

P-i  F 


“ 

nl 

v„ 

— B 

= 

-B 

Jl 

Hi 

- 

. 

....(102) 


We  now  wish  to  examine  the  properties  of  the  matrices  and  in 

particular  ask  whether  U,  is  symmetric  and  whether  any  are  related  to 
■p  n  K  P 

each  other.  Now  =  Z  ^pkZp  an<*  a^ter  some  rearrangement  we  obtain 


rj  ^ 

u.  =  y  s ,  z . 

k  L  pk  p 
p=i  v  v 


....(103) 
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if  m 


2q  are  symmetric;  the  rest  are 


Hence,  only  Uj  and  U^+J 
unsymmetric.  Next  let  us  consider 


U 


T 

n-r+ 


2 


n 


K 


L  p,n-r+2 

p=i 


Z  . 
P 


After  some  manipulation  we  find 


U 


T 

n-r+2 


=  U 

r 


(104) 


..(105) 


So  for  storage  purposes  we  only  require  space  for  the  matrices  (1/n)  Z_, 

D 

Z.,  0lt  U  2,  ...  U  where  r  is  the  integer  part  of  n/2.  Of  these  matrices 
(l/n)  Zg,  Uj  and,  for  n  =2q,U^+j  are  symmetric.  Thus,  no  more  storage  is 
required  for  this  set  of  matrices  than  for  the  original  set  and  the  two 
sets  can  share  the  same  storage  area  since  the  original  set  are  never 
required  again  once  the  U.  are  calculated. 


All  these  transformed  matrices  are  factorised  either  by  MA29 

if  they  are  symmetric  or  by  MA23  if  they  are  not.  A  new  routine,  MA23T, 

T 

has  been  written  to  solve  for  the  system  A  x  =  ^  where  the  factorization 
of  the  unsymmetric  matrix  A  is  known,  since  the  Harwell  package  does  not 
include  this  possibility. 


The  only  other  problem  concerns  the  triangle  functions  which 
overlap  from  one  section  to  the  next.  Such  triangle  functions  must  exist 
otherwise  each  section  will  be  electrically  isolated  from  the  other.  The 
geometry  package  in  the  CHAOS  program  has  been  modified  to  provide  these 
overlaps  automatically.  To  this  end  a  convention  has  been  adopted  in 
which  these  overlaps  are  attached  to  the  left-hand  side  of  a  section. 
This  is  illustrated  in  figure  10(a).  A  further  modification  to  the  CHAOS 
geometry  package  is  required  to  ensure  that  wires  attached  to  the  axis  of 
symmetry  have  their  overlap  connections  arranged  symmetrically.  This  is 
done  by  ensuring  that  each  such  wire  has  an  overlap  on  to  the  axis,  as 
il'ustrated  in  figure  10(b).  This  means  that  if  a  set  of  wires  meet  on 
the  axis,  then  there  must  be  another  wire  lying  along  the  axis  on  to 
which  the  set  of  wires  must  overlap.  If  such  an  axis  wire  does  not  exist, 
then  due  to  the  method  of  attaching  overlaps  (described  in  section  4.1), 


a  symmetric  system  of  overlaps  cannot  be  set  up.  This  problem  could  be 

got  round  by  placing  the  unsymmetric  set  of  triangle  functions,  due  to 

this  junction  on  the  axis,  into  the  matrix  Z  .  However,  this  has  not  been 

□ 

programmed,  as  it  would  cause  further  problems  in  that  the  column  of 
matrices  Z^  in  matrix  (82)  would  have  to  be  replaced  by  a  column  of 
unequal  matrices.  This  would  wreck  the  analysis  that  has  been  described 
here . 


7 . 2  Left- right  symmetry 

This  section  is  included  although,  at  the  time  of  writing,  the 
option  is  not  part  of  the  CHAOS  program.  A  plane  must  be  chosen  to  be  the 
mirror  plane.  As  it  is  envisaged  that  aircraft  will  be  one  of  the  more 
important  bodies  for  which  this  facility  is  useful,  it  is  assumed  that 
the  mirror  plane  is  the  XZ  plane.  This  does  not  affect  the  analysis  given 
here  but  will  affect  the  computer  programming. 


The  matrix  takes  the  form 


(106) 


where  the  square  matrix  Z  now  represents  the  interactions  between 

D 

triangle  functions  lying  wholly  in  the  mirror  plane  or  which 
symmetrically  straddle  the  plane.  The  square  matrix  Zl  represents  the 
interactions  between  triangle  functions  on  one  side  of  the  mirror  plane, 
while  the  square  matrix  Z2  represents  the  interactions  between  triangle 
functions  on  one  side  of  the  plane  with  those  on  the  other  side.  Lastly, 
the  rectangular  matrix  Z represents  the  interaction  between  triangle 
functions  in  the  plane  (or  straddling  the  plane)  and  those  on  one  side  of 
the  plane.  Once  again  the  analysis  is  restricted  to  symmetric  impedance 
matrices  so  that  storage  is  saved  and  the  principle  of  reciprocity  is 
satisfied.  The  manner  in  which  equation  (106)  is  written  implies  that  the 
numbering  scheme  on  one  side  of  the  boundary  is  a  reflection  of  that  on 
the  other.  This  implies  too  that  the  matrix  Z2  is  itself  symmetric,  so 
that 
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T 

Z2  "  Z2 


....(107) 


The  corresponding  current  and  voltage  vectors  are  written  as 


'i  ' 

V„  " 

-B 

-B 

I  = 

Ii 

and  V  * 

lx 

I2 

1* 

....(108) 


We  define 


J,  =  Ii 

+  I2, 

Wx  =  Vj  +  V2 

J2  =  Ij 

"  12, 

W2  -  Vj  -  v2 

Ux  =  zx 

+  z2 , 

U2  =  Zx  -  Z2 

...(109) 


then  the  equations  we  wish  to  solve,  analogous  to  equations  (100)  and 
(102),  are 


U2J2  -  W2 , 


....(110) 


1  Z 

2  B 

r 

« 

»— *1 

CM 

t _ 

r 

>1 

i _ 

ZA 

a 

1 _ 

lx 

— i 

dTi 

_ i 

We  then  obtain 

11  -  \  ( Ji  +  J2), 

12  =  I  ( Jx  -  J2 ) . 


....(Ill) 


....(112) 


Now  since  Z\  and  Z2  are  symmetric,  Ui  and  U2  are  symmetric  and  advantage 
should  be  taken  of  this  to  save  storage.  As  the  matrix  in  equation  (111) 
is  symmetric,  the  Harwell  routine  MA29  may  be  used  for  all  the  factori¬ 
zations. 
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I  Two  planes  of  mirror  symmetry 

Once  again  this  section  is  provided  although  the  programming 
has  not  yet  been  done.  It  is  included  because,  for  an  aircraft  sitting  on 
a  perfectly  conducting  ground  plane,  we  have  both  the  left-right  symmetry 
of  the  aircraft  and  its  reflection  in  the  ground  plane.  This  method  would 
be  invalid  for  a  ground  plane  of  finite  conductivity.  Following  section 
7.2,  the  left- right  plane  may  be  thought  of  as  the  XZ  plane  and  the  ground 
plane  as  the  XY  plane,  but  this  does  not  affect  the  analysis  given  below. 


The  matrix  takes  the  form 


Z  = 


T 

ZA 

*1 

zT 

zT 

Zi 

zj 

zj 

z  I 

z2 

Zi 

T 

Z4 

zj 

m 

z3 

Z4 

Zi 

7T 

z2 

Z2 

z3 

z2 

Zi . 

understood 

by  reference  to  figure  11.  Yet 

....(113) 


square  matrix  Z  represents  the  interactions  between  triangle  functions 
lying  wholly  in  one  or  other  mirror  plane  or  which  symmetrically  straddle 
one  or  other  plane.  The  square  matrix  Zj  represents  the  interaction 
between  triangle  functions  lying  exclusively  in  one  of  the  four 
sections.  The  square  matrices  Z2,  Z3  and  Z4  represent  the  interaction 
between  triangle  functions  lying  wholly  in  section  1  with  those  lying 
wholly  in  sections  2,  3  and  4  respectively.  They  also  represent 
permutations  of  the  above,  eg,  Z2  also  represents  the  interaction  between 
sections  3  and  4,  Z3  between  2  and  4,  and  Z*  between  2  and  3.  Once  again 
the  numbering  within  each  section  must  follow  the  numbering  in  the  first 
section.  In  this  way  the  four  matrices  Zi,  Z2,  Z3  and  Z4  are  all 
themselves  symmetric  provided  that  the  original  matrix  Z  is  symmetric. 
Lastly,  the  rectangular  matrix  Z^  represents  the  interactions  between 
triangle  functions  in  the  mirror  planes  (or  straddling  it)  and  one  of  the 
four  sections. 
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In  view  of  the  above  we  have 

Z£  ■  Z.f  i  ■  1,  2,  3,  4.  ....(114) 


The  corresponding  current  and  voltage  vectors  are  written  as 


ri  * 

rv_ 

-B 

-B 

II 

Vi 

I2 

and  V  “ 

V2 

13 

V3 

1 4 

V4 

In  order  to  solve  for  the  currents  we  define 


Ji  *  Ii  +  I2  +  I3  +  I*  ) 

“  "  "  "  } 

J2  =  Ij  -  I2  +  I3  -  I4  ) 

“  “  “  } 

J3  -  +  I2  ~  I3  -  1%  } 

"  "  } 

J4  =  -  j[2  -  +  I4  ) 


ana  Wx  =  Vx  +  V2  +  V3  +  V„  } 

"  “  “  } 

Wz  -  Vi  -  V2  +  V3  -  V*  } 

”  } 

w3  =  Vi  +  V2  -  V3  -  V4  } 
"  "  } 
W4  =  Vi  -  v2  -  V3  +  } 

and  also  Ui  =  Zx  +  Z2  +  Z3  +  ZH  } 

} 

U2  =  Zi  ~  Z2  +  Z3  ~  Zi*  } 

} 

U3  =  Zj  +  Z2  “  Z3  ~  Z4  } 

} 

Uij  ■  Zj  -  Z2  “  Z3  +  Z%  • } 


We  now  solve  the  equivalent  equations 


U.J.  -  W.,  i  -  2,  3,  4 

1—1  — 1  ’  ’  ' 


(115) 


....(116) 


....(117) 


•  see 


(118) 


(119) 


42 


-£  z 

4  B 

*:  ■ 

1 

-p* 

.  ZA 

Ui 

£1 

L*  J 

and 


Lastly,  we  obtain  the  currents  from  the  transformations 


£i  "  T  (£1  +  J2  +  J3  +  Ji*)  } 

1  "  “  } 

12  =  r  (Ji  -  J2  +  J3  ~  J4)  } 

1  "  ) 

13  =  T  (Jl  +  J2  -  J3  -  J4)  } 

7  } 

£<♦  =  r  (£1  ~  £2  ~  £3  +  £4 )  •  } 


....(120) 


....(121) 


Once  again  we  note  that  Ult  U2 ,  U3  and  U4  are  all  symmetric  as  is  the 
matrix  in  equation  (120).  This  can  be  taken  advantage  of  to  reduce 
storage  requirements  and  allows  the  Harwell  routine  MA29  to  be  used  for 
all  the  factorizations. 


The  major  programming  difficulty  in  both  sections  7.2  and  7.3 
will  come  in  the  geometry  package,  in  setting  up  the  correct  overlaps  for 
wires  crossing  symmetry  planes  and  in  adding  in  the  correct  contributions 
to  finite  conductivity  -  this  latter  is  probably  the  most  difficult  but 
it  is  anticipated  that  much  of  the  programming  provided  for  the 
rotational  symmetry  option  in  section  7.1  will  carry  over  in  a 
straightforward  manner.  If  two  wires  meet  at  the  symmetry  plane,  then  the 
overlap  will  need  to  be  put  in  the  matrix  Z^,  unless  there  are  three 
segments  under  the  triangle  function  on  one  side  of  the  symmetry  plane 
and  one  segment  on  the  other.  When  more  than  two  wires  meet,  then  to 
obtain  a  symmetric  set  of  overlaps,  it  will  be  necessary  for  one  of  these 
wires  to  lie  in  the  symmetry  plane. 

8.  SPARSE  MATRIX  APPROXIMATION 


Inevitably  there  are  problems  which  require  matrices  too  large 
to  fit  into  the  computer.  One  scheme  for  reducing  the  storage  needed  for 
bodies  a  few  wavelengths  long  is  to  approximate  the  full  matrix  by  a 
sparse  matrix  by  neglecting  interactions  between  parts  of  the  body  which 
are  more  than  a  certain  distance  apart. 
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A  special  sparse  matrix  version  of  CHAOS  has  been  written 
incorporating  this  idea.  It  uses  the  Harwell  routine  MA28  [13],  modified 
for  complex  arithmetic,  to  factorize  the  impedance  matrix,  Z,  into  the 
sparse  factors  which  are  upper,  U,  and  lower,  L,  triangular.  Thus,  we 
have 


Z  -  L.U. 


(122) 


Unfortunately,  there  is  no  reliable  routine  available  to  produce  a 
symmetric  factorization  of  Z  although  Harwell  are  expecting  to  produce 
one  in  due  course  (it  is  provisionally  named  MA27).  This  has  the 
consequence  that  more  storage  is  used  than  is  necessary.  Also  it  has  been 
noticed  that  in  the  L  and  U  factors,  there  are  many  small  elements  which 
can  be  neglected  without  significantly  affecting  the  radar  cross- 
sections.  Now,  because  the  factorization  is  unsymmetric,  the  dropping  of 
these  small  elements  breaks  the  symmetry  and  so  the  principle  of 
reciprocity  is  not  exactly  satisfied.  So  there  are  at  least  two 
compelling  reasons  for  changing  to  a  symmetric  factorization  routine  as 
soon  as  possible. 


9.  COMPLEX  FREQUENCY  VERSION  OF  CHAOS 


Let  us  consider  the  equation  (29)  which  may  be  written 
symbolically  as 


ZI  =  V, 


....(123) 


where  Z  is  the  known  impedance  matrix,  V  is  the  known  voltage  vector  and 
1^  is  the  unknown  current  vector.  Now  the  complex  frequencies,  S^,  at 
which  the  determinant  of  Z  is  zero  are  given  by 


det  Z(S.)  -  0. 
i 


•  •  •  • 


(124) 


The  current,  for  some  complex  frequency  S,  may  be  expanded  as 


I(S) 


(125) 
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In  the  time  domain,  we  have  the  Fourier  transform 
“  S.t 

I(t)  =  l  e  1  R..  _ (126) 

i=i  -1 

So  the  establish  the  resonance  frequencies  and  the  decay  times  of  the 
structure  and  depend  on  the  structure  only.  They  are  independent  of 
polarisation  and  orientation  of  the  structure.  However,  the  amplitudes 
are  dependent  on  polarisation,  orientation  and  structure,  but  not  on 
frequency.  This  scheme  is  known  as  the  singularity  expansion  method  114] 
and  is  abbreviated  to  SEM  in  the  literature. 

Each  structure  is  thought  to  have  a  unique  set  of  S.  but  it 
requires  a  considerable  amount  of  computing  time  to  calculate  each  S.. 
The  R^  may  be  calculated  rather  more  easily  once  the  S.  are  known.  In 
certain  circumstances,  the  can  be  extracted  from  experimental 

measurements  of  the  time  dependent  response  of  a  structure  to  a  Gaussian 
or  wide  band  pulse. 

The  CHAOS  program  has  been  modified  to  calculate  the  complex 
resonant  frequencies  S...  It  uses  a  search  procedure  in  the  complex 
frequency  plane  due  to  Muller  115]  .  This  method  essentially  takes  three 
guesses  to  a  zero  and  fits  a  parabola  through  the  values  of  the 
determinant  of  the  corresponding  impedance  matrices.  The  next  approxi¬ 
mation  to  the  zero  is  found  as  a  zero  of  this  parabola.  The  first  guess  is 
then  discarded  and  the  approximation  just  found  is  substituted.  The 
process  is  then  repeated  until  convergence.  However,  this  scheme  is  not  a 
very  efficient  one  as  it  can  fail  to  find  some  zeroes  and  has 
difficulties  with  double  zeroes.  Its  main  advantage  is  that  it  does  not 
require  the  derivative  of  the  impedance  matrix  to  be  calculated. 

How  the  initial  three  starting  values  are  chosen  needs  some 
explanation.  We  need  to  study  equations  (125)  and  (126)  and  remember 
that  we  are  interested  in  real  solutions  to  real  problems.  This  implies 
that  the  occur  in  complex  conjugate  pairs  and  must  lie  in  the  left- 
hand  half  of  the  complex  frequency  plane  in  order  to  produce 
exponentially  decaying  currents.  For  convenience,  we  choose  to  work  with 
normalised  frequencies  and  define 
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....(127) 


where  c  is  the  velocity  of  light  and  L  is  the  maximum  linear  dimension  of 
the  scattering  object.  With  these  units  our  starting  guesses  for  a  zero 
are:- 


S(l)  -  -  0.1  +  i0.5  +  t, 

=  -  0.15  +  t,  _ (128) 

^(s)  -  -  0.1  +  t, 

where  t  may  be  supplied  as  data  to  the  program  if  there  is  any  good  reason 
for  knowing  what  the  zero  is.  If  it  is  not  supplied  as  data,  it  is  taken 
to  be  zero.  In  the  case  of  a  dipole  one  might  take  t  to  be  i0.5,  il.O, 
il.5,  etc,  corresponding  to  the  known  resonances  of  the  dipole.  These 
values  have  been  found  by  trial  and  error  to  be  sensible.  There  is  one 
final  problem  that  needs  mentioning;  the  impedance  matrix  Z(S)  is 
singular  when  the  frequency  is  zero.  Thus,  in  the  computer  program,  one 
works  with  the  combination  SZ(S)  which  is  finite  at  the  origin. 

A  better  scheme  based  on  contour  integration  has  been  devised 
by  Baum  et  al.  (16  ]  .  The  idea  is  to  divide  the  complex  frequency  plane  up 
into  boxes  and  to  evaluate  a  certain  integral  around  the  edge  of  this 
box.  A  set  of  equations  is  obtained,  one  for  each  integral,  which  define 
how  many  zeroes  lie  inside  the  box  and  provide  a  set  of  simultaneous  non¬ 
linear  equations  to  solve  for  the  required  complex  frequencies.  If  there 
are  too  many  zeroes  in  the  box,  it  is  a  trivial  task  to  further  subdivide 
it.  The  equations  that  are  being  used  are 

M 

2?l^csk  Is  log  Z(S)dS  *  l  si>  k  “  0,  1  ...  M,  ....(129) 

where  M  is  the  number  of  zeroes  inside  the  contour  C.  This  scheme  has  not 
been  programmed  into  CHAOS  but  should  these  zeroes  be  required  in  future, 
it  would  be  sensible  to  change  over  to  Baum's  method.  It  is  unlikely  to 
reduce  the  computing  time  very  much. 
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FIGURE  1 


STEP  FUNCTION  APPROXIMATION  TO  TRIANGLE  FUNCTION  Tj 
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FIGURE  4 
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(B)  DEGENERATE  CASE 


WIRE 


FIGURE  6  Position  of  the  Field  Point 

co 


WIRE  1 


(a) 


OVERLAPPING  WIRE  JUNCTION  TREATMENT  USED  BY  CHAO  AND  STRAIT. 
(OVERLAPPING  TRIANGLE  FUNCTIONS  SHOWN  DASHED.) 


WIRE  1 


(bi  JUNCTION  TREATMENT  WITHOUT  OVERLAPS 


FIGURE  7 
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-  WIRE  OVERLAP 

.  OUTLINE  OF  TRIANGLE  FUNCTION 

(  i  )  WIRE  NUMBER  i 

FIGURE  8  Triangle  Function  i  Leaves  Junction  J. 

Straddles  Junction  K  and  Enters  Junction  L 


Case  (i)  a^  (1)  ”  0 
b.  (1)  =  1 

i 


Case  (ii)  a.  (1)  ■  0,  (1)  “  0, 

fa.  (2)  *  0,  a.  (3)  -  0,  b.  ,  (3)  =  l] 
1  1-1  1-1  1-1  J 


l  '  '  ’  l 

b.  (2)  =  1 

i 


[a._j  (4)  =  0,  b._i  (4)  =  0] 


Case  (iv)  (2)  =  0,  a^  (3)  -  0, 
b.  (3)  =  1 

[a.+i  (1)1=  0,  b.+l  (1)  »  0] 


TRIANGLE  FUNCTION  i 
ENTERS  JUNCTION  L 


Case  (v) 

a.  (4)  =*  0,  b.  (4)  =  0 

K+1  <2>  ■  «i+l  (3>  -  bi+1  «>  * 


b.  (4)  -  1 


FIGURE  9  Various  Situations  at  a  Junction 
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■  WIRE 
•  OVERLAP 


FIGURE  10(a)  Illustration  of  Convention  Adopted  for  Adding  Overlaps  to  a  Section 
to  Cause  the  Whole  Structure  to  be  Electrically  Connected 


FIGURE  10(b)  Illustration  of  Arrangements  of  Overlaps 
Connecting  Wires  to  Axis  in  a  Symmetric  Fashion 
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T 


FIGURE  11 


Section  Numbering  Scheme  for  2  Planes  of  Mirror  Symmetry 
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APPENDIX  A 


ACCURACY 


A1 . 


EVALUATION  OF  FUNCTION 


Let  us  write  the  function  of  equation  (49)  as 


.  1  rA/2  e 

*  =  4TfS  U/2  “ R—  dZ  * 


-ikR 


.  (A1 ) 


where  we  have  chosen  a  local  cylindrical  coordinate  system  based  on  the 
source  segment  as  Z  axis.  The  segment  length  is  taken  as  A  and  the  origin 
of  coordinates  is  at  the  centre  of  the  segment.  Here  R  is  given  by 


R  =  /p2  +  (Z  -  Z')2, 


. . (A2) 


where  p  and  Z  are  the  cylindrical  coordinates  of  the  field  point  in  this 
system.  The  geometry  is  illustrated  in  figure  A1 .  Since  we  have 
stipulated  that  the  field  point  lies  on  the  wire  surface,  while  the  line 
source  lies  on  the  wire  axis,  the  integrand  can  never  become  singular. 


For  r  <  5A,  Harrington  recommends  the  expansion 
*  »  |l,  -  ik(  I  j  -  rl,)  -  y-  (I,  -  2rl2  +  r*I,)  ♦  i  ^  (I*  -  3rl,  +  3r*I2  -  r3I,)j, 


....(A3) 


where  Ij  =  iog 
I2  =  A, 


Z  + 


A/2  +  /p2  +  (Z  4-  A/2) 


n 


l-z 


-  A/2  +  /p2  +  (Z  -  A/2)s 


.... (A4) 

.... (A5) 


Z  till  /p2  +  (Z  +  A/2)2  +  i2llf-Z2  /p2  +  (Z  -  A/2)2  +  §—  Ix , 

. . . . (a6) 


I3  =  - ^ 1 —  /p-  T  t  - - 2 - r  w  x"  '  2 


I  =  Ap2  +  3A3 /I 2  +  AZ2 


i . (A7) 
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While  for  r  *  5  A,  he  recommends 


II 

|rt> 

d  ? 

i - * 

> 

o 

+ 

i£ 

> 

+ 

(£)\ 

♦  i(^)3Aa 

*  (“)'*.] , 

00 

< 

• 

• 

• 

where 

A  = 
o 

■  *  |  (f;)f  i  *  3(f)2] 

■r[3  -  3o( 

:§)’  *  35(f)'], 

.... (A9) 

Ax  = 

i  &[- 1  *  3(f)!] 

-  30(-)2 
r 

'  35(f)']  , 

....(A10) 

A2  = 

-  i  -  n  - 

12(f)2  * 

15(f)'], 

....(All) 

A3  = 

w  C§r)[3C|52  -  »(f)']. 

.... (A12) 

Ah  * 

T5o  <?)*• 

.... (A13) 

A2 .  MAXIMUM  SEGMENT  LENGTH 


We  next  need  to  know  how  accurate  these  formulae  are.  Both 
equations  (A3)  and  (A8)  are  the  result  of  integrating  power  series 
expansions  term  by  term  and  truncating  after  the  fourth  and  fifth  terms 
respectively.  The  ratio  of  the  fourth  term  to  the  leading  term  in 
equation  (A3)  is  (n/N)3  and  the  ratio  of  the  fifth  term  to  the  leading 
term  in  equation  (A8)  is  ( tt/N)  h  where 


N 


_  X  _  2tt_ 
~  A  kA’ 


(A14) 


ie,  N  is  the  ratio  of  wavelength  to  segment  length.  The  use  of  higher 
order  terms  in  these  expansions  has  been  investigated  and,  although  the 
accuracy  of  the  impedance  matrix  elements  is  improved,  the  accuracy  of 
the  radar  cross-sections  is  less  affected  for  a  given  value  of  N.  The 
accuracy  of  the  matrix  elements  for  various  values  of  N  is  indicated  in 
table  A1 . 
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TABLE  A1 


Percentage  Accuracy  of  Matrix  Elements 


Formula  for  ip  Function 

(A3) 

(A8) 

10 

3 

1 

7 

9 

4 

6 

14 

8 

5 

25 

16 

4 

48 

38 

In  CHAOS,  the  ratio  N  is  printed  out  for  each  wire  with  a 
warning  if  it  is  less  than  10.  If  it  is  less  than  6  for  any  wire,  the  job 
is  rejected.  Just  because  the  matrix  elements  are  accurate,  it  does  not 
mean  that  the  current  distributions  and  radar  cross-sections  are 
accurate.  For  simple  geometries,  like  a  dipole,  one  may  need  as  many  as 
50  segments  per  wavelength  especially  if  one  is  near  a  resonance  and 
wishes  to  get  cross-sections  accurate  to  a  few  per  cent.  This  is  because 
the  frequency  at  which  a  dipole  resonates,  and  the  cross-section  at  this 
frequency,  change  by  about  10%  as  the  number  of  segments  is  increased 
from,  say,  10  segments  per  wavelength  to  50.  Increasing  this  figure 
further  leads  to  only  small  changes  of  the  order  of  1%.  However,  because 
the  resonance  is  so  narrow,  if  one  is  interested  in  a  spot  frequency 
close  to  the  resonance,  the  cross-section  may  be  wrong  by  100%  if  one 
uses  10  segments  per  wavelength.  For  wires  made  of  material  of  low 
conductivity  (say,  of  the  order  of  10"  mho/m),  for  which  the  resonance 
effect  is  not  very  pronounced,  fewer  segments  are  necessary. 

For  complicated  geometries,  it  appears  that  one  can  get  away 
with  14  segments  per  wire  for  good  conductors  and  10  segments  per  wire 
for  other  materials.  This  should  give  the  positions  of  lobes  in  the  radar 
cross-section  pattern  correct  to  about  ±  5°  and  the  amplitudes  of  the 
largest  lobes  correct  to  about  ±  2  db.  It  is  not  possible  to  be  more 
precise  because  of  the  impossibility  of  using  shorter  segments  due  to 


computer  storage  limitations  and  because  more  accurate  experimental  data 
on  complicated  geometries  does  not  appear  to  be  available. 


It  is  up  to  the  user  to  satisfy  himself  that  his  calculations 
are  properly  converged  by  staging  calculations  with  more  segments  where 
possible,  checking  against  whatever  experimental  data  he  may  have  and 
using  his  common  sense. 

A3.  SMALL  SEGMENT  LENGTHS 


Intuitively,  one  would  expect  to  get  a  better  and  better  answer 
to  a  calculation  as  the  segment  length  is  reduced.  This  is  not  found  to 
be  the  case  and  needs  explanation.  Let  us  consider  a  straight  dipole  of 
length  L  and  evaluate  the  matrix  elements  given  by  equation  (32)  in  the 
limit  as  A  Che  segment  length  tends  to  zero  keeping  the  wavelength  and 
radius  fixed.  We  find  after  some  manipulation  that 


Z. . 
Ji 


where 


Z  A  + 
0) 


-ikr 


■nrr 

P 


(3  +  3ikr  -  k2r2) 
_ P.  P 

iuer2 


(10  if  j  =  i, 

( 

-  (  3  if  j  -  i  ±  1, 

J  ( 

(  0  otherwise, 

p  =  2(j  -  i), 


r  =  /a2  +  PW. 
P 


•  •  •  a 


(A15) 


The  matrix  order  for  such  a  problem  is  (L/2  A)  -  1.  The  structure  of  this 
impedance  matrix  is  curious.  The  main  diagonals  and  the  diagonals  on 
either  side  of  it  are  or  the  order  of  A  for  very  small  segment  lengths 
while  all  other  terms  are  of  the  order  of  A2.  In  any  one  row  of  this 
matrix  there  are  roughly  L/2  A  elements  of  the  order  of  A 2  so  that  these 
could  add  up  to  give  a  terra  of  the  order  of  A  which  would  then  be 
comparable  to  the  diagonal  term.  It  has  been  found  that  when  a  32  bit 
computer  word  length  is  used  problems  start  to  arise  if  there  are  1000  or 
more  segments  on  a  wire.  Thus,  the  CHAOS  program  will  reject  all  problems 
which  contain  any  wire  divided  into  more  than  500  segments. 
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Let  ub  now  recall  what  is  happening  physically  to  our  segment 
as  this  limit  is  approached.  It  starts  out  being  in  the  form  of  a  long 
thin  cylinder  and  changes  through  a  short  fat  cylinder  to  eventually 
become  a  wafer  thin  disc.  The  formulae  used  in  CHAOS  are  based  on  the 
assumption  of  a  long  thin  cylinder  and  we  need  to  know  how  short  the 
cylinder  can  become  before  this  assumption  is  violated.  Let  us  define  the 
aspect  ratio  R  to  be  the  ratio  of  segment  length  to  radius.  From 
numerical  experiments  with  CHAOS  using  a  dipole,  we  find  that  when  R  is 
as  small  as  0.1,  the  current  distributions  become  highly  oscillatory  with 
spatial  period  2A.  These  current  distributions  appear  to  be  all  right  for  R 
equal  to  0.5.  Thus,  we  have  decided  to  reject  all  jobs  for  which  R  is  less 
than  0.5  but  warning  messages  are  printed  out  if  R  is  less  than  5. 

A4.  ANOTHER  LIMITING  CASE 


Another  limiting  case  which  is  of  interest  is  that  which  occurs 
when  the  segment  length  A  tends  to  zero  in  such  a  way  that  the  ratio  of 
segment  length  to  radius,  R,  remains  constant  at  the  same  large  value. 
After  some  further  manipulation  we  find  the  impedance  matrix  elements  to 
be 


Zjj  4irajeA 


2  l0g  a  ~  V 


*  r  2  A  +  0(kA) , 
4  u) 


Zj,j±i  4TTweA 

z  =  i{ie"iqA 
j ,  i  AircoeA 


K  AQ  \ 

*{v 


.... (A1 6) 

.... (Al 7 ) 


2(q2  -  1)  q2  -  4  2(q2  -  9) 


+  ikA  { - - -  +  — £ -  +  - 1 - } 

2(q2  -  1)  q2  -  4  2 ( q 2  -  9) 

where  i^j.  j  ±  1  and  q  *  2  |j  -  i  | . 


+  0(kA),  . . (A18) 


Suppose  now  we  restrict  our  dipole  to  have  four  segments  and  to 
be  a  perfect  conductor,  then  the  radar  cross-section  will  be 


a0f>m 


Kk^L6  sin**  9 


. . . . (A19) 


T 


ve)  ■  °’ 


(9)  -  0 


where  L  =  4  A  and  0  is  the  angle  between  the  direction  of  the 
transmitter/receiver  and  the  axis  of  the  dipole.  <j)  is  the  azimuthal  angle 
about  this  axis.  The  notation  CL,  means  the  radar  cross-section  for  the 

Dip 

transmitter  polarised  in  the  direction  of  increasing  0  ,  while  the 
receiver  is  polarised  in  the  direction  of  increasing  <p- 


From  the  book  by  Ruck  [17],  the  radar  cross-section  for  a 
perfectly  conducting  short  thin  cylinder  of  length  L  and  radius  a  is 


V9)  ■ 


TTk-L6  sin4  6 
144  (log  -  l]2' 


. . . .(A20) 


V0)  ■  °- 


0.^(6)  =  0. 


These  two  sets  of  formulae  are  not  quite  the  same  and  reflect  the  fact 
that  in  CHAOS  the  current  distribution  is  forced  to  be  triangular  in 
shape  while  the  formulae  given  by  Ruck  are  obtained  by  treating  the 
dipole  as  the  limiting  case  of  a  prolate  spheroid.  These  are  only  equal 
when  L/a  is  about  157  and  for  larger  values  of  L/a  the  CHAOS  prediction 
will  be  smaller  than  those  of  Ruck. 


A5.  PROBLEMS  AT  VERY  LOW  FREQUENCIES 

The  behaviour  of  the  impedance  matrix  as  the  frequency  tends  to 
zero  is  yet  another  interesting  limit.  In  this  case  we  can  expand  e 
in  equation  (45)  as 

e_lkR  =  1  -  ikR  +  0(k2 ) ,  - (A21) 


since  kR  is  small  compared  to  1  when  the  frequency  is  so  low  that  the 
largest  dimension  of  the  scattering  body  is  much  smaller  than  the 
wavelength.  Thus, 


ZE.  = 


ji  4TTiu)e 


i—  f  dl  f  dd-  LL.  T.(d)  <  T.U') 


wires  wires 


R  dd  j  dd'  i' 


where  the  term  of  order  unity  integrates  to  zero. 


+  0(k) , 
.... (A22) 
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From  equation  (44)  we  also  have 


=  0(k) , 


. .(A23) 


while  from  equation  (47)  using  equations  (78)  and  (79)  we  have 


0(1). 


(A24) 


Thus,  the  imaginary  part  of  the  impedance  matrix  is  0(1 /k)  while  the 


real  part  is  0(1)  for  non-perfect  conductors  and  0(k)  for  perfect 
conductors.  So  we  can  split  the  matrix  into  three  parts  as 


Z..  =  R..  +  iL..  +  T.., 
Ji  Ji  Ji  Ji 


where  R. . 

Ji 

L. . 
Ji 

T.  . 
J1 


We  wish  to  solve 


*  0(k), 

=  o(I), 

=  0(1). 


. . . . (A25) 


(R  +  iL  +  T) I  =  V, 


.... ( A26 ) 


where  =  0(1)  as  k  ■+■  0,  for  the  currents  1^.  In  general  this  should 
present  no  problems  and  at  low  frequency  is  equivalent  to  solving 


iLI  =  V 


. . . . (A27) 


to  give  currents  of  0(k)  which  is  intuitively  sensible  since  at  low 
frequency,  the  whole  structure  is  at  the  same  potential  and  no  currents 
will  flow.  The  trouble  comes  for  those  structures  for  which  the 
determinant  of  L  is  singular,  as  in  the  case  of  a  planar  loop  of  wire.  For 
this  case  equation  (A27)  has  many  solutions  and  for  a  given  finite 
computer  word  length  there  is  a  frequency  below  which  the  problem  cannot 


Thus,  a  test  has  been  programmed  into  CHAOS  to  detect  the 

situation  where  the  matrix  Z  is  nearly  singular.  This  is  done  by 

T 

examining  the  factorized  form  of  Z,  be  it  LDL  (for  symmetric  matrices) 
or  LU  (for  unsymmetric  cases).  The  maximum  and  minimum  elements  of  the 
diagonals  of  either  D  or  II  are  found  and  the  ratio  is  taken.  From  tests  on 
a  dipoLe,  it  has  been  found  that  when  this  ratio  is  greater  than  25000 
for  a  computer  word  length  of  32  bits,  the  results  will  be  spurious.  In 
CHAOS,  if  the  ratio  exceeds  20000  the  calculation  is  terminated  and  no 
results  can  be  obtained.  In  cases  where  this  ratio  is  only  just  less  than 
20000  (the  maximum  and  minimum  are  printed  out  by  the  program),  results 
should  be  checked  to  see  that  they  are  intuitively  sensible.  For 
computers  which  use  a  64  bit  word  length,  the  results  appear  to  be 
sensible  for  values  of  this  ratio  up  to  100000;  beyond  this  value  they 
progressively  deteriorate. 

A6.  COMPUTER  WORD  LENGTH  PROBLEMS 


The  CHAOS  program  has  been  made  to  work  on  both  an  IBM  computer 
using  32  bit  word  length  and  a  CRAY  computer  using  64  bit  word  length. 
Tests  have  been  carried  out  on  both  machines  using  a  matrix  of  order  650. 
Tn  general,  no  sigificant  change  in  results  has  been  noticed.  In 
particular,  at  and  near  the  peaks  of  a  radar  lobar  diagram,  there  is 
usually  no  change  to  four  significant  figures  although  occasionally  a 
change  of  one  unit  in  the  fourth  figure  does  occur.  In  the  troughs 
between  the  lobes,  the  changes  can  be,  and  are,  much  bigger  and  radar 
cross-sections  which  are  six  orders  of  magnitude  smaller  than  the  maximum 
using  32  bit  word  length  can  become  fourteen  orders  of  magnitude  smaller 
using  a  64  bit  word  length.  However,  these  are  large  changes  in  small  and 
insignificant  numbers  and  so  are  of  no  consequence. 

The  only  other  change  is  that,  as  already  mentioned  at  the  end 
of  the  previous  section,  the  use  of  a  longer  word  length  allows 
calculations  to  be  done  at  slightly  lower  frequencies  than  can  be  done 
with  a  32  bit  word  length. 
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FIGURE  A1  Geometry  for  Colculat 


APPENDIX  B 


CONVENTIONS  USED  IN  CHAOS 


BI.  INTRODUCTION 


This  appendix  details  the  time  and  phase  conventions,  together 
with  the  coordinate  systems  that  may  be  used  to  describe  the  effect  of  an 
electromagnetic  wave  being  scattered  by  a  body.  The  system  as  described 
here  is  used  in  the  CHAOS  computer  program. 

B2 .  WHY  USE  A  COMPLEX  REPRESENTATION  FOR  A  REAL  EFFECT? 


The  first  point  to  explain  is  why,  when  we  are  dealing  with 
real  electric  fields  and  real  currents  in  a  real  world,  we  should 
•represent  these  by  complex  quantities.  The  short  answer  is  that  this  is 
mathematically  convenient  as  may  be  seen  by  studying  the  structure  of 
Maxwell's  equations  for  the  electromagnetic  fields. 

From  equations  (1),  (3)  and  (58),  it  may  be  seen  that 

V-V-6  +  ua  fr  £  +  ue  ~  £  *  0,  - (Bl) 

3t  -  3t2  - 

where  z ,  u,  and  a  are  assumed  to  be  constant  and  S  is  the  real  electric 
field.  Equation  (Bl)  can  be  further  simplified  by  assuming  a 
cosinusoidal  time  variation  at  a  constant  frequency  U).  We  could  write  all 
quantities  in  the  form 

^.(r,t)  =  Aj(r)  cos  (cut  +  40  for  j  =  x,y,z,  ....(B2) 

but  substitution  into  equation  (Bl)  then  gives  an  inelegant  and  rather 
cumbersome  equation  to  solve.  It  is  much  more  convenient  mathematically 
to  write 


E(r,t)  =  E(£)e*Wt  > 


. . (B3) 
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where  E^(£)  =A^(r)e  for  j  =  x,y,z,  . ,..(B4) 
and  then  to  solve  V,*V~E(£)  -  T2E(£)  *  0,  . ...(B5) 
where  T2  =  -  iuBi(o  +  icse).  . ,..(B6) 

Once  we  have  the  solution  to  the  complex  equation  (B5)  it  is  a 
trivial  job  to  obtain  the  solution  (B2)  to  the  real  problem. 


B3.  TIME  CONVENTION  AND  COORDINATE  SYSTEM 


Following  the  discussion  of  section  B2  the  time  variation  will 
be  described  by  e1U)t,  where  co  is  the  frequency  in  radians/second.  The 
complex  time  dependent  electric  field  is  defined  as 

E(r,t)  =  E(£)e1(i,t.  - (B7) 

The  real  time  dependent  electric  field  is  given  by 


|(£»t)  =  ReE(£,t).  _ (B8) 

When  a  radar  looks  at  an  object,  there  are  in  general  two 
coordinate  systems  which  are  of  interest.  The  first  is  a  system  fixed 
with  respect  to  the  object  as  used  in  the  computer  program  CHAOS.  In  this 
system  the  Z  axis,  of  a  right-handed  cartesian  system  OXYZ,  might  for 
instance  be  chosen  to  be  the  axis  about  which  the  body  has  some  natural 
symmetry.  In  this  system  the  radar  beam  would  be  incident  from  a 
direction  defined  by  ( 0,  <f> ) ,  the  usual  spherical  coordinates  (see 
figure  Bl(a)).  Thus,  the  incident  electric  field  transverse  to  the 
direction  of  propagation  would  have  components  in  the  9  and  <J> 
directions. 

The  second  system  of  coordinates  is  fixed  with  respect  to  the  radar 
and  is  again  a  right-handed  cartesian  coordinate  system  Oxyz  in  which  the  radar 
wave  propagates  along  the  z  axis  from  large  positive  z  towards  z  *  0.  It  is 
convenient  to  take  the  x  axis  in  the  direction  of  increasing  9,  and  y  in  the 
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direction  of  increasing  <J»,  where  these  directions  are  as  defined  in  the 
first  system.  In  figure  Bl,  these  axes  are  offset  along  the  radar  line  of 
sight  for  clarity.  In  the  program  CHAOS,  it  is  assumed  that  the  object 
axis  (OZ)  is  vertical  so  that  the  y  axis  is  necessarily  horizontal.  The  x 
axis  then  forms  the  third  member  of  an  orthogonal  triad.  Conventionally 
(as  in  the  nomenclature  "VV  polarisation")  the  x  axis  is  referred  to  as 
"vertical"  though  this  is  a  misnomer  unless  the  radar  line  of  sight  lies 
in  the  horizontal  plane.  For  this  special  case  (figure  Bl(b))  Ox  is  truly 
vertical,  actually  vertically  downwards.  It  would  be  more  sensible  if 
CHAOS  print  outs  used  the  phrase  8  and  4>  polarisation  in  place  of 
vertical  and  horizontal  polarisation. 

In  the  radar  coordinate  system,  a  typical  complex  electric 
wave,  travelling  inwards  from  large  z  towards  z  =  0  along  the  z  axis, 
would  be  represented  by 

i(kz+<}>  )  i(kz-Hj>  ) 

E(r)  =  (A  e  X  ,  A  e  ^  ,  0),  ....(B9) 

—  —  x  y 

where  k  =  2tt/X  is  the  wave  number,  X  =  2ttc/d  is  the  wavelength  and  c  is 
the  velocity  of  light.  The  sign  of  z  is  changed  for  a  wave  travelling 


to  the  time  origin  and  only  the  phase  difference  is  of  any  physical 
significance.  The  two  amplitudes  A  and  A  could  be  replaced  by  an 

•  /~ 2  r  .  x  y, 

amplitude  y  A^  +  A^  times  the  appropriate  cosine  or  sine  of  some  angle; 
however,  it  is  simpler  to  treat  them  in  this  way. 


B4 .  POLARISATION  CONVENTIONS 


The  real  time  dependent  electric  field  is 
i(u)t+kz+<$>  )  i(wt+kz+<|>  ) 

£(r,t)  =  Re(Axe  *  ,  A^e  y  ,  0).  ....(BIO) 


This  vector  lies  in  the  xy  plane  and  is  in  general  rotating  about  the  z 
axis.  The  angle  that  £  makes  with  the  x  axis  is  ¥  where 


tan  ¥ 


A  cos  (u)t  +  kz  +  A  ) 

y _  y 

A  cos  (oj t  +  kz  +  d>  )’ 
X  Tx 


....(Bll) 
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The  magnitude  of  &  is 


l£l  =  /a*  cos2  (tot  +  kz  +  <b  )  +  A2  cos2  (u>t  +  kz  +  6  ).  - (B12) 

A  X  y  '  y 

For  linear  polarisation'}'  is,  by  definition,  constant.  Since 


, .  A  A  sin  (<b  -  <b  ) 

dij*  _  x  y _ ^x_*_y 


. . . . (B1 3) 


in  general  this  implies  <{>-<))=  mr,  where  n  is  any  integer.  The  electric 

x  y 

vector  therefore  makes  a  constant  angle  ±  tan  (A  /A  )  with  the  x  axis. 

y  x 

The  other  two  possible  solutions  to  equation  (B13)  are  A  =0  (vertical 
polarisation)  and  A^  =  0  (horizontal  polarisation)  and  these  are  also 
just  special  cases  of  the  general  solution. 


If  0  <  <}>  and  A  and  A  have  the  same  sign,  then  4'  decreases 

x  y  x  y 

as  time  increases  and  elliptic  polarisation  in  a  clockwise  (right-hand) 
direction  about  the  direction  of  propagation  and  looking  along  the 
direction  of  propagation  is  obtained*  (see  figure  B2).  The  polarisation 
is  elliptic  because  the  magnitude  of  the  electric  field  given  by  equation 
(B12)  traces  out  an  ellipse  in  the  xy  plane. 

In  order  to  obtain  circular  polarisation  the  magnitude  of  8 
must  be  independent  of  time.  For  this  to  occur,  A  must  equal  A  and  |>  = 

x  y  y 

^x  +  ^2n  *  Dir/ 2  where  n  is  any  integer. 


It  is  also  instructive  to  ask  what  shape  the  beam  has  at  some 


fixed  instant  of  time. 


For  t  constant, 


dY 

d(kz) 


A  sin  (d>  -  d)  ) 

y _  x _ y  _ 


A  cos2  (wt  +  kz  +  d>  ) 

X  X 


(B14) 


If  *  <  <b  and  A  and  A  have  same  sign,  then  ¥  decreases  as  one  moves 
x  y  x  y  ° 

along  the  wave  (from  large  z  to  small  z)  in  the  direction  of  propagation, 
so  that  a  left-hand  "corkscrew"  is  obtained.  Thus,  this  left-handed 
(anti-clockwise)  "corkscrew"  in  space  moves  forward  and  rotates  in  a 
right-handed  (clockwise)  direction  as  time  increases  when  looking  in  the 
direction  of  propagation. 


*This  is  the  right-hand  thread  rule  referred  to  the  direction  of 
propagation  which  agrees  with  the  definitions  of  the  IRE  1942  standards  [18] 


B5.  INCIDENT  AND  SCATTERED  FIELDS 


The  incident  field  will  be  assumed  to  be  a  plane  wave  of  the 
type  given  by  equation  (B8).  This  will  be  correct  provided  the 
transmitter  is  sufficiently  far  from  the  scattering  obstacle. 

Using  the  asymptotic  form  for  the  scattered  field,  it  can  be 
shown  that  the  field  transverse  to  the  direction  of  propagation  is  0(l/r) 
where  r  is  the  distance  of  the  field  point  from  the  scatterer  while  the 
field  in  the  direction  of  propagation  is  0(l/r2).  So  at  large  distances 
from  the  scatterer  the  field  becomes  essentially  transverse  to  its 
direction  of  propagation.  The  behaviour  of  the  components  of  the  field 
determine  the  polarisation  of  the  scattered  wave.  The  polarisation  of  the 
receiver  will  determine  which  of  these  components,  if  any,  is  measured. 

B6 .  SCATTERING  MATRIX 


An  incident  elliptically  polarised  field  EX(r)  can  be 
expressed  in  terms  of  cartesian  coordinates  (as  shown  in  section  B3)  and 
so  may  be  related  to  cartesian  coordinates  of  the  backscattered  field 

g 

E  (r)  by  a  scattering  matrix  119],  S,  where 

ES(r)  =  SEX(r).  _ (B15) 

Having  found  the  scattering  matrix,  S,  for  a  given  pair  of 
orthogonally  polarised  fields,  it  is  convenient  to  be  able  to  transform 
this  matrix  to  S1,  where  S1  refers  to  a  different  pair  of  orthogonally 
polarised  fields  (eg,  S  might  refer  to  linearly  polarised  fields  in 
the  0  and  <p  directions  and  S1  to  circularly  polarised  fields  with  left- 
and  right-handed  rotations  about  the  direction  of  propagation).  This  may 
be  done  by  means  of  a  change  of  polarisation  basis  and  the  fields  may  be 
expressed  in  terms  of  an  elliptic  basis  instead  of  a  linear  one  by  a 
coordinate  transformation  U. 

Now  equation  (B9)  may  be  re-written  as 

“  Vx  +  E^y>  . . . .  (B16) 
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where  a  and  u  are  unit  vectors  in  the  x  (or  9)  and  y  (or  <$>)  directions. 

—x  — y 

Suppose  u^  and  u^  are  an  arbitrary  pair  of  orthogonal  vectors  coplanar 

with  u  and  u  ,  then  they  must  be  obtained  from  them  by  a  unitary 
—x  y 

transformation.  In  this  case  the  matrix  U  has  the  property  that 

U*TU  =  I,  (B17) 

where  T  denotes  transpose  and  *  denotes  complex  conjugate. 

So  the  incident  electric  field  in  the  elliptic  basis,  is 

given  by 

I1  =  u£l(r).  _ ( B18) 

Then  ES(r)  =  SU_1UE1(r)  =  SU-1^.  _ (B19) 

The  scattered  field  is  transformed  by  the  same  transformation 

so  that 

Z.S  =  UES(r)  =  USlfV,  (B20) 

s 

where  is  the  scattered  vector  in  the  elliptic  basis.  This  gives  the 
scattering  matrix  S1  =  USU  1  referred  to  an  elliptic  basis. 

Suppose  the  transformation,  U,  in  equation  ( B 1 8 )  gives  a 
right-hand  elliptically  polarised  incident  electric  field,  then  the 
transformation  given  by  equation  (B20)  produces  a  left-hand  elliptically 
polarised  scattered  field  since  the  direction  of  propagation  is 
reversed.  This  is  because  when  we  talk  of  left-  or  right-handedness  a 
direction  is  associated  with  it  and  when  we  speak  in  this  way  we  are 
effectively  using  a  different  cartesian  system  for  the  two  directions  of 
propagation  -  one  being  right-handed  for  the  incident  wave,  the  other 
left-handed  for  the  scattered  wave  as  viewed  by  an  observer  at  large 
distance  (eg,  at  the  radar). 

It  is  convenient  to  use  a  different  transformation  on  the 
scattered  field  so  that  in  the  example  above  a  right-handed  elliptically 
polarised  scattered  field  is  produced.  We  can  do  this  by  using 


=  U*ES, 


. .. . (B21 ) 


so  that  the  scattering  matrix 

S1  =  U*SU_1  _ (B22) 

gives  right-hand  transmit,  right-hand  receive  and  left-hand  transmit, 

left-hand  receive  on  the  main  diagonals  of  S1  .  By  taking  the  complex 

conjugate  of  U,  we  have  effectively  changed  the  signs  of  $  and  $  in 

x  y 

equation  ( B 1 1 ) ,  thus  changing  the  sense  of  the  polarisation. 

B7.  FORM  OF  TRANSFORMATION  MATRIX 


The  unitary  condition  (B17)  has  been  imposed  on  U  to  preserve 
magnitudes  of  vectors,  so  only  four  of  the  eight  quantities  (the  four 
complex  u ^ j )  in  the  transformation  are  independent.  The  most  general  such 
matrix  may  be  written 


mn 

cos  ^e  , 


mi  2 


-  sin  £e 


i(^i\+n22“Hv2)  it>22 

sin  £e  ,  cos  £e 


. .(B23) 


The  determinant  of  U  1  is  e*  ^  1+r|2  ^  .  Thus,  the  determinant  of 
S  is  not  preserved  under  the  transformation  (B22)  but  has  its  phase 
altered.  We  follow  other  authors  [20-22]  in  imposing  the  condition  that 
the  determinant  of  the  transformation,  U,  is  unity.  The  most  general  such 
matrix  is  now 


if1 


inn  mi2 

cos  £e  ,  -  sin  £e 

-ini2  -inn 

sin  £e  ,  cos  £e 


.  (B24) 


This  contains  three  unknowns  which  may  be  related  to  the  phase, 
ellipticity  and  the  orientation  of  the  elliptical  vector  formed  by  the 
transformation  [21]  .  This  may  be  brought  out  by  making  the  following 
substitutions:- 
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tan  rii!  =  -  tan  tan  a, 
tan  n 1 2  =  ~  cot  4>  tan  a> 

„  1  +  tan2  \b  tan2  a 

tan  x,  =  tan-1  \|»  - - - 1 — ~  • 

tan2  4>  +  tan  a 

The  transformation  matrix  [22]  now  becomes 

Q  =  U_1  =  R(4>)H(a),  - (B26) 


where  R(40 


and  H(a) 


),  -  sin  (J/l 
>,  cos  4)  ] 


cos  - 
sin 

cos  a,  i  sm  a" 
i  sin  a,  cos  a 


. . (B27) 


. . . . (B28) 


We  shall  go  on  to  show  that  R(<JO  is  a  rotation  matrix  that  determines  the 
orientation  of  the  major  axis  of  the  ellipse  and  that  H(  a  )  is  an 
ellipticity  matrix  which  determines  the  ratio  of  major  to  minor  axes.  We 
nave  reduced  the  3  unknowns  in  equation  (B24)  to  2  in  equation  (B26) 
effectively  by  choosing  an  arbitrary  phase  factor  but  this  does  not 
matter  since  phase  can  only  be  determined  to  within  an  arbitrary  additive 
constant.  The  extra  condition  we  have  imposed  can  be  written  as 


tan  £ 


cos  n12 
cos  Tit  I 


/ tan  rti  i  cot  ru  2  » 


(B29) 


but  this  appears  to  have  no  physical  significance. 


It  is  worth  noting  that  the  matrices  R  and  H  do  not  commute, 
ie,  RH  ^  HR.  With  the  transformation  matrix  given  by  equation  (B26)  we 
see  that  the  determinant  of  S  is  preserved  under  the  transformation  [22] 

*t 

since  Q  Q  =  I  and  det  Q  =  1  so  that 

det  S1  =  det  (QTSQ)  =  det  QT  det  S  det  Q  =  det  S.  _ (B30) 


Also  preserved  is  the  trace  of  the  power  scattering  matrix  [2J.  This 
matrix  is  defined  as 


P  =  S*TS. 


• (B31 ) 
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Now 


pi  = 


1  *T  1 

S  S  , 


*T  *T  *  T 

=  (Q  s  XQ  ) (Q  SQ), 

*T  *T 

=  Q  S  SQ, 

*T 

=  Q  PQ. 


...(B32) 


The  trace  (ie,  the  sum  of  the  diagonal  elements)  represents  the 
total  power  that  would  be  returned  to  a  pair  of  orthogonally  polarised 
antennae  and  is  given  by 


Now 


1  S1 1 

1 2  +  1  S22 1 

2  +  1 S12 ( ' 

=  p 

(Q*  Q  + 

Q*  Q  )  + 

1 1 

11  11 

12  12 

+  p 

(Q*  Q  + 

Q*  Q  )  + 

2  1 

2  1  11 

22  12 

=  P 

+  P  , 

1  1 

22 

II 

H 

^0 

N-K» 

.... (B33) 


12  11  21 

>  (Q*  Q 

22  21  21 


12  22 

)*  q  : 
22  22 


where  the  P. .  and  Q. .  are  the  elements  of  P  and  Q  respectively. 

ij  ij 


B8. 


ROTATION  MATRIX 


For  positive  the  matrix  R(\p)  operating  on  the  vector  V  will 
rotate  V  in  the  xy  plane  by  an  amount  'p  in  a  counter  clockwise  (left- 
hand)  direction  looking  along  the  direction  of  propagation  (see  figure 
B3).  Following  a  rotation  the  elements  of  the  scattering  matrix  become:- 


*1  1 

=  S,  ! 

cos2 

<P 

+  (s 

1 2 

♦  S2 1 ) 

COS  Ip 

sin 

<p  + 

S22 

sin2 

't'j 

S 1  2 

=  (S22 

-  S 

ll) 

cos 

+ 

sin  \J; 

+  S12 

2 

cos4 

*  - 

S2  1 

sin2 

sji 

=  (S22 

-  s 

11) 

cos 

* 

sin  ip 

+  S2i 

2 

COS 

- 

Si  2 

sin2 

'P, 

S22 

=  S1 1 

sin2 

* 

~  (S 

1 2 

*  S2i) 

COS 

sin 

4*  + 

S22 

cos2 

Note  that  these  do  not  agree  with  reference  [22]  in  which  a  sign  has  been 
changed  between  their  equation  (12)  and  their  equation  (18). 
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We  see  that 


Su  +  S22  J 
Si  2  -  S2 1  . 

B9 .  ELLIPTICITY  MATRIX 


and 


si  1 


+  s 


l 

22  " 


si  2 


si  i 


The  matrix  H(  a  )  operating  on  the  vector  (  q  )  produces  the 

vector  (jCg£n  ^  )•  Using  equation  (B9)  this  may  be  interpreted  as  an 

electric  field  for  which  A  =  cos  a,  A  «  sin  a  and  <{>  =  0,  <t>  =  tt/2.  The 

x  y  x  y 

magnitude  of  the  field  is  (from  equation  (B12)) 

ISI  s2  a  cos  (ait  +  kz)  +  sin2  a  sin2  (ait  +  kz), 

so  that  2  traces  out  an  ellipse  with  ratio  of  minor  to  major  axes  given  by 
tan  a.  Provided  a  is  positive,  then  the  circulation,  looking  in  the 
direction  of  propagation,  is  clockwise  (or  right-handed)  as  shown  in 
figure  B4(a)*.  When  a  =  ±  Tt/4,  we  get  circular  polarisation. 

When  H(a)  operates  on  (  °  )  an  elliptically  polarised  field  is 
again  produced,  the  major  axis  now  being  at  right-angles  to  that  obtained 
from  H(a)(  q  )  and  the  circulation  direction  being  of  opposite  sense,  as 
shown  in  figure  B4(b).  Thus,  the  phase  zero  of  H(a)(  J  )  leads  that  of 
H(a)(  J  )  by  90°. 


Suppose  H(a)  operates  on 

so  that  A  =  A  =  — - ,  =  <fe  =  a. This  is  a  linearly  polarised 

x  y  /2  x  y 

field  (figure  B4(c)).  This  should  not  be  surprising  since  we  have  noticed 
that  the  rotation  and  ellipticity  matrices  do  not  commute  and  to  obtain 

r  1  1 

[/2  j 

elliptic  polarisation  we  are  expecting  if  we  operate  with  H(  a  )  on  either 
(  J)  or  (  1  )>  ie,  before  operating  with  the  rotation  matrix. 

♦Note  that  this  does  not  agree  with  reference  [22]  unless  their  conventions 
on  circulation  direction  are  opposite  to  ours. 


/2 

1 


from  (J)  requires  a  rotation  of  45°.  Thus,  we  shall  only  get  the 


is 


1_ 

ia 

✓2 

e 

1__ 

ia 

/2 

e 

1_ 

/2 

1_ 

1/2  J 


then  the  vector  field  obtained 
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Following  a  change  of  ellipticity  the  scattering  matrix 
elements  become :- 


s{, 

=  S(  1  cos2  a  + 

i  ( si  2 

+  S2i) 

sin  a  cos 

a 

- 

S22 

sin2 

a, 

s}2 

=  i ( S 1 1  +  S2  2  ) 

sin  a 

cos  a  + 

S12  cos2 

a 

- 

S2  1 

sin2 

a, 

S2  1 

=  i ( Si  1  +  S2  2  ) 

sin  a 

cos  a  - 

S12  sin 

a 

+ 

S2 1 

cos2 

a, 

S22 

=  Six  sin2  a  + 

i(S12 

+  S2 1 ) 

cos  a  sin 

a 

+ 

S22 

cos2 

a. 

Note  that  1  +  Sj2 

=  ( S4 1 

+  S22 ) 

cos  2a  + 

i ( Si  2  + 

S21 ) 

sin 

si  1  -  s2 2 

=  Stl 

-  S22  > 

sj2  -  S2 1 

=  Si  2 

-  S2 1  » 

....(B36) 


and  that  the  phase  zero  of  s}2,  S21  leads  sjj  by  90°  while  that  of  S22 
leads  S } !  by  180°. 

For  circular  polarisation  a  =  +ir/4  for  right-hand  and  a  =  -tt/4 
for  left-hand  rotation.  This  gives  (fora  =  +tt/4):- 


sii(n/4)  =  -  Sj2  +  +  S2\))» 

S*  2  (tt/4)  =  j  (  i  ( S 1 1  +  S22)  +  Si2  ~  S21), 

Sz  1  (tt/4)  =  (  i  ( S 1 1  +  S22)  “  S12  +  S2i)» 

si2(Ti/4)  =  4  (~  Su  +  S22  +  i(S12  +  S21)). 


..(B37) 


We  may  change  the  notation  slightly  to  make  this  more 

meaningful.  By  remembering  our  conventions,  we  have  for  right-hand 

transmit  and  right-hand  receive  S  i  1  (tt/4)  which  we  re-write  as  C_  where 

RR 

the  first  suffix  refers  to  the  receiver  and  the  second  suffix  to  the 
transmitter.  Similarly  for  transmit  in  the  direction  of  increasing  9  and 
receive  in  the  direction  of  increasing  0  we  have  Si  1  which  we  re-write  as 
Sqq.  In  this  way  equation  (B37)  may  be  rerwritten  as:- 
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T 


Crr  "  2  (  S<H>  +  S90  +  l(StO  +  S0<J>) }  ’ 


crl  ~  2  (i(s<h>  +  S^0  +  se<t> 5 * 

clr  =  I  (i(s<t>4>  +  W  +  s<i>9  "  W’ 


. . . . (B38) 


C  =  ~  (+  s..  -  Saa  +  i(S.0  +  SQ , ) ) , 
LL  2  'p<p  00  q>0  0<p  ’ 


where  again  the  first  suffix  refers  to  the  receiver  and  the  second  to  the 

transmitter.  The  phase  zero  of  C  and  C  leads  that  of  C  by  90°  while 

KL  LR  RR 

that  of  C  leads  that  of  C  by  180°. 

L  L  RR 


BIO.  COMBINED  EFFECT  OF  ROTATION  AND  ELLIPTICITY  CHANGE 


Since  the  rotation  matrix  and  the  ellipticity  matrix  do  not 
commute,  the  transformation  (B35)  must  be  performed  first  followed  by 
equation  (B36).  The  expressions  for  the  elements  of  the  new  scattering 
matrix  in  terms  of  the  original  are  not  particularly  informative  and  so 
are  not  given  explicitly.  However,  for  the  special  case  of  circular 
polarisation  we  have:- 


C1  =  c  e"2i^ 
RR  RR  ’ 


C1  =  C  e21*, 
LL  LL  ’ 


CLR  CLR’ 


. . . .(B39) 


C 


RL* 


Despite  the  differences  noted  earlier,  ref.rence  [22]  agrees  with  these 
formulae . 

We  now  define  four  radar  cross-sect' ons  invariant  under 
rotations  and  these  are:- 


'’8 


■  ■  1  I  CB0  !  2  *  I VI'  *  lCR,|2  *  lC,J2J. 


j  -  W  I CLR| 2 , 


K  ■  4,t>  IC^I1, 


L  -  4,r»  ICJ", 


. . . . ( B40) 


assuming  the  incident  plane  wave  has  unit  modulus.  Here  r  is,  as  before, 
the  distance  from  the  scatterer  to  the  receiver.  Note  that 


1  <|seel!  *  lsJ!  *  lv|!  *  |se*|! 


. .. .(B41) 


also,  and  that  for  a  monostatic  radar  (ie,  one  for  which  the  transmitter 
and  the  receiver  are  in  the  same  place) 


I  =  2J  +  K  +  L. 


. . . .(B42) 


For  a  body  of  revolution,  with  the  coordinate  axis  along  the  body  axis, 
we  have  in  the  monostatic  case  [23] 


S0<t>  =  °’  S*0  =  ° 


....(*43) 


and  for  arbitrary  rotation  \p  of  the  polarisation  axes:- 


1 

(s 

\  - iV 
-  S  Je 

2 

00  <t>4> 

1 

(S  -  s  )e2i*, 

2 

00  <M> 

=r 1  = 

lrl 

:  1  i(S00  +  W’ 

(B44) 


SQ0  =  SO0  cos2  *  *  sin2  t, 

S0<t>  =  S<£0  =  (S<J><j>  '  S90)  cos  Sln 
=  S00  sin2  *  +  S4,4>  cos2  *• 


,(B45) 
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Bl  l. 


DEPOLARISATION 


Bickel  [22]  defines  a  quantity  D  ,  which  he  calls  the 


depolarisation,  as 


D  =  1 

n 


S  v  S.  12 

'if) 


2Tr(P 


2  |  S96  Sjxt>'  *  I  S94> | 2  *  |  S(j)9  1  2 
^  s9o  1 2  +  lse<tJ2  +  ^  s4>e  ^ 2  +  ls04)l2 


,  (B46) 


'ClJ2  +  iCRR!2 


|c  I2  +  |c  I2  +  I c  I2  +  |c  I2 

[  LL1  1  LR1  1  RL1  !  RR1 


A  sphere,  for  which  C  and  C  are  zero,  has  D  =0  while  a 

LtLj  kk  ^  n 

dipole,  for  which  Sg^  =  S^g  =  =  0,  say,  has  Bickel  then  goes 

on  to  show  that  larger  values  of  characterize  bodies  with  multiple 
scattering  centres.  is  invariant  to  rotations  but  not  to  ellipticity 
changes.  Dn  is  calculated  by  CHAOS  automatically  but  its  usefulness  has 
not  yet  been  assessed. 
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APPENDIX  C 


LIST  OF  SYMBOLS 


Many  of  these  symbols,  when  used  in  the  text,  have  suffices  or 
superfices  associated  with  them.  Not  all  of  these  are  separately  defined 
here.  This  list  is  not  complete  but  covers  the  most  frequently  used 
symbols. 


Cl.  COORDINATE  SYSTEMS  AND  GEOMETRY 


OXYZ 

Oxyz 

r.  =  (x,y,z) 
t 

a 

£ 


0,<t> 


— X 


-y 

R  =  |r-r' | 


ds 

ri 

d£ 

T(£) 

Din 


A 


object  coordinate  system 
radar  coordinate  system 
general  coordinate  vector 
time 

wire  radius 
length  along  wire 

spherical  polar  coordinates  in  OXYZ  system 
unit  vectors  in  Ox  and  Oy  directions 

distance  between  observation  point,  £,  and  some  source 
point,  jr' 

element  of  surface  area 
unit  normal  to  surface 
element  of  wire  length 
triangle  function 
unit  vector 

scalar  product  between  vectors  in  directions 

r.  .  -  r.  and  r.  .  -  r.  where  r.  is  the  vector  from 
— jm+1  — jm  —  m+1  —in  — jm 

the  origin  to  the  start  of  segment  m  under  triangle 

function  T. 

~3 

segment  length 


C2. 


ELECTRIC  FIELDS 


g(r,t) 

E(r,t) 

E(r) 

D 

H 

B 

J 

1(4) 

P 

o(4) 

e 


y 

a 

-  _i_ 

v 


u)  =  2ttv 


k  =  2-rr/A  =  - 
c 


T  =  /-iwy (a+iwe) 
A 

4> 

A. 

J 

♦j 

¥ 

e 


real  electric  field 
complex  electric  field 

complex  electric  field  assuming  e1Wt  time  variation 

complex  electric  flux 

complex  magnetic  field 

complex  magnetic  flux 

complex  conduction  current 

complex  line  density  of  current 
complex  electric  charge 
complex  line  density  of  charge 
dielectric  constant 

magnetic  permeability 

conductivity 

velocity  of  light  in  medium 

frequency  in  cycles/sec 

frequency  in  radians/sec 

impedance  of  medium 

wavelength 

wave  number 

complex  wave  number 

magnetic  vector  potential 

electric  scalar  potential 

magnitude  °f£j>  j  =  x»y>z 

phase  of  j  =  x,y,z 

angle  between^  (x,t)  and  the  x  axis 

complex  electric  field  intensity  referred  to  an  elliptic 
polarisation  basis 
voltage  vector 
impedance  of  wire 


z 

g(0.« 

p 

<p 

s 

n 

R,L,C 


*  impedance  matrix 

*  radar  cross-section  in  direction  (0,<J>) 

*  power 

=  integral  defined  by  equation  (49)  or  (Al) 
=  skin  depth 
=  loss  angle 

*  resistance,  inductance  and  capacitance 


C3.  SCATTERING  MATRICES 


S 

*T 

P  =  S  S 
U 


?,nll,nl2,ri22 

a 

i> 

R(H>) 

H(a) 

Q  =  U-1  =  R(i|»)H(a) , 

see,se<}>’s<i>9’s<M> 


CR R’  CRL’  CLR’  °LL 


scattering  matrix,  elements  S.^ 
power  scattering  matrix,  elements 
unitary  transformation  matrix 

the  four  independent  quantities  in  the  unitary  matrix  U 
ellipticity  angle  (arctan  of  ratio  of  minor  to  major 
axes  of  ellipse) 

rotation  angle  (of  body  about  radar  line  of  sight) 

rotation  matrix 

ellipticity  matrix 

elements  Q. . 

ij 

are  the  elements  of  the  scattering  matrix  S  referred  to 
linearly  polarised  fields  in  the  directions  of  increasing 
0  or  <f>.  The  first  suffix  referring  to  the  receiver,  the 
second  to  the  transmitter 

are  the  elements  of  the  scattering  matrix  S  referred  to 
circularly  polarised  fields  with  right,  R,  or  left,  L, 
hand  circulation 
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47rr2  {  |  CRR 1  2  +  |  C 


+  C. 


♦  K.i.r 


!clrI2 

!crrI2 


47rr 


^LL 1 


RL1  1  LR1  1  LL1 
the  sphere  accept  radar  cross-section 

a  sphere  reject  radar  cross-section 

another  sphere  reject  radar  cross-section 

depolarise  ion. 
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So—  Metric  and  SI  Unit  Conversion  Factor* 

(Baaed  on  DEF  STAN  00-11/2  "Metric  Units  for  Uae  by  the  Mlnletrv  of  Defence", 
DS  Met  5501  "AWRE  Metric  Guide"  and  other  British  Standards) 


Quantity  Unit  Syabol  Conversion 


Basic  Units 


Length 

metre 

m 

1  m  -  3.2808  ft 

1  ft  -  0.3048  m 

Mass 

kilogram 

kg 

1  kg  -  2.2048  lb 

1  lb  -  0.45359237  kg 

1  ton  -  1016.05  kg 

Derived  Units 

Force 

newton 

N  •  kg  m/s2 

1  N  -  0.2248  lbf 

1  lbf  -  4.44822  N 

Work,  Energy,  Quantity  of  Heat 

joule 

J  -  N  m 

1  J  -  0.737562  ft  lbf 

1  J  -  9.47817  «  10‘4  Btu 

1  J  -  2.38846  »  10~4  kcal 

1  ft  lbf  -  1.35582  J 

1  Btu  -  1055.06  J 

1  kcal  -  4186.8  J 

Power 

watt 

U  -  J/e 

1  W  -  0.238846  cal/s 

1  cal/s  -  4.1868  W 

Electric  Charge 

coulomb 

C  •  A  a 

- 

Electric  Potential 

volt 

V  •  W/A  -  J/C 

- 

Electrical  Capacitance 

farad 

F  -  A  s/V  -  C/V 

- 

Electric  Resistance 

ohm 

n  -  v/a 

- 

Conductance 

8  lenten 

S  -  1  !)"' 

- 

Magnetic  Flux 

weber 

Wb  •  V  b 

- 

Magnetic  Flux  Density 

tesla 

T  -  Wb/m2 

- 

Inductance 

henry 

H  -  V  s/A  -  Wh/A 

Complex  Derived  Units 

Angular  Velocity 

radian  per  second 

rad/s 

1  rad/s  •  0.159155  rev/s 

1  rev/s  »  6.28319  rad/s 

Acceleration 

metre  per  square  second 

m/s2 

1  m/s2  -  3.28084  ft/s2 

1  ft/s2  -  0.3048  m/s2 

Angular  Acceleration 

radian  per  square  second 

rad/s2 

- 

Pressure 

newton  per  square  metre 

N/ra2  -  Pa 

1  N/m2  -  145.038  *  10-&  lbf/ii 
1  lbf /In?  -  6.89476  *  10*  N/m: 

bar 

bar  -  10 5  N/m2 

1  in.  Hg  -  3386.39  N/m2 

Torque 

newton  metre 

N  m 

1  N  m  -  0.737562  lbf  ft 

1  lbf  ft  -  1.35582  N  m 

Surface  Tension 

newton  per  metre 

N/m 

1  N/m  -  0.0685  lbf /f t 

1  lbf /f t  -  14.5939  N/m 

Dynamic  Viscosity 

newton  second  per  square  metre 

N  s/m2 

1  N  s/m2  -  0.0208854  lbf  s/ft 
1  lbf  s/ft2  -  47.8803  N  s/m2 

Kinematic  Viscosity 

square  metre  per  second 

m2/s 

1  m2/s  -  10.763°  ft2/s 

1  f t 2 / s  •  n.0929  m2/s 

Thermal  Conductivity 

watt  per  metre  kelvin 

W/m  K 

- 

Odd  Units* 

Radioactivity 

becquerel 

Bq 

1  Bq  -  2.7027  «  10~1!  Cl 

1  Cl  -  3.700  «  1010  Bq 

Absorbed  Dose 

gray 

cy 

1  Gy  -  100  rad 

1  rad  -  0.01  Gy 

Dose  F.qul valent 

a lever t 

Sv 

1  Sv  •  100  rem 

1  rem  -  0.01  Sv 

Exposure 

coulomb  per  kilogram 

C/kg 

1  C/kg  -  3876  R 

1  R  -  2.58  «  lO"4  C/kg 

Rate  of  leak  (Vacuum  Systems) 

millibar  litre  per  second 

mb  1/s 

I  mb  *  0.750062  torr 

1  torr  -  1.33322  mb 

*These  teres  are  recognised  terns  within  the  netrlc  svatea. 


